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We analyze how a short distance boundary condition for the Schrodinger equation must change as 
a function of the boundary radius by imposing the physical requirement of phase shift independence 
on the boundary condition. The resulting equation can be interpreted as a variable phase equation 
of a complementary boundary value problem. We discuss the corresponding infrared fixed points 
and the perturbative expansion around them generating a short distance modified effective range 
theory. We also discuss ultraviolet fixed points, limit cycles and attractors with a given fractality 
which take place for singular attractive potentials at the origin. The scaling behaviour of scattering 
observables can analytically be determined and is studied with some emphasis on the low energy 
nucleon-nucleon interaction via singular pion exchange potentials. The generalization to coupled 
channels is also studied. 
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I. INTRODUCTION 

One of the most remarkable features of low energy scat- 
tering for a short range spherically symmetric potential 
is the onset of scale invariance and rotational invariance. 
For a short range potential, U (r) , which will be the sub- 
ject of the present paper, with a typical size a and in the 
long wavelength limit fca <C 1, the Schrodinger's equa- 
tion for the reduced wave function in the I — channel 
becomes 



u"{r)^0. 

The solutions to this equation are straight lines 



u{r) 



br . 



(1) 



(2) 



where the coefficients a and b (or rather their ratio), 
are determined by matching to the solution in the re- 
gion where the potential acts. The intersection of the 
asymptotic solution with the x-axis determines the scat- 
tering length a = —a/b. Eq. ([1]) is obviously invariant 
under the scale transformation r Xr, but the general 
solution is not due to the presence of the potential U{r). 
There are, however, two particular situations where the 
solution also transforms well under scaling, namely when 
either a = or 6 = 0. These two cases correspond to 
a = (trivial scattering) and a = ±00 (zero energy 
bound state) respectively. The presence of a potential 
U{r) at r ^ a or finite energy, fca ~ 1, induces scaling 
violations which can be computed within perturbation 
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theory, and obey scaling laws. As a matter of fact, it is 
interesting to see what kind of perturbation theory can be 
constructed around these two very simple scale invariant 
cases. To answer this question Wilsonian renormalization 
group (RG) methods seem the adequate tool 

Our interest and focus in the present paper is mainly 
motivated by encouraging developments in the last 
decade in nuclear physics, and more specifically on the 
Nucleon-Nucleon (NN) interaction problem in the frame- 
work of effective field theories triggered by Weinberg's 
work [1, 0| (for reviews see e.g. Refs [3, M, S S [I])- 
Actually, much of the discussion has been unavoidably 
linked to the ability of designing adequate regularization 
schemes which in addition to preserve the symmetries can 
be removed beyond perturbation theory, allow to han- 
dle highly singular potentials and provide a one-valued 
renormalization group flow of low energy parameters. 

In this paper we attack the problem by using the 
boundary condition (BC) regularization. It is based 
on the natural idea that all unknown information be- 
low some scale R can always be parameterized in terms 
of a mixed boundary condition at the distance R (see 
Sect.|n|. Strictly speaking, this is exactly true for a non- 
relativistic system which can be described through the 
Schrodinger equation, since it is a second order differen- 
tial operator. The long distance physics is assumed to 
be determined in terms of a potential U (r) above the 
boundary radius R. The BC method has been used ex- 
tensively in the past for the treatment of NN scatter- 
ing [9j|, and more recently for analyzing the commonly 
accepted analysis of phase-shifts in terms of the chiral 
expansion 11| but always keeping the short distance cut- 
off finite. The equivalence between effective field theory 
(EFT) and boundary conditions at the origin for short 
range potentials has been established in Ref. We as- 
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sume this equivalence to hold also in the presence of long 
distance interactions. Actually, we will see that it is pos- 
sible to shrink the boundary to the origin with a smooth 
limit in the physical observables. From a Lagrangian 
viewpoint non-local momentum dependent terms in the 
interaction can, after suitable field redefinitions and us- 
ing the equations of motion, be expressed as local and 
center-of-mass (CM) energy dependent potentials 

The boundary condition satisfies a renormalization 
group equation (see Sect. IIII[) which infrared and ultra- 
violet critical points may be examined and character- 
ized. Moreover, the boundary condition may be inter- 
preted in terms of an outer truncated potential problem 
(see Sect. IIV[) . The renormalization group equation for 
the boundary condition becomes a variable phase equa- 
tion p^ . which has a straightforward physical interpre- 
tation. Besides, one of the virtues of the variable phase 
approach is that it always deals with a on-shcU problem, 
i.e., at any stage of the calculation the variable phase shift 
exactly corresponds to a physical phase shift of a certain 
on-shell problem. The price to pay is the non- linear char- 
acter of the equation. Thus, the well known off-shell am- 
biguities characteristic of Lippmann-Schwinger equation 
(LSE), and which make the discussion on renormalization 
cumbersome, particularly when truncations are involved, 
never appear. 

Although some of the aspects presented in this pa- 
per have been known or implicitly assumed for a long 
time in some way or another, we believe that the inter- 
pretation in terms of a short distance boundary condi- 
tion becomes quite transparent and unifying, particularly 
when explicit long range, i.e. non-contact, interactions 
are considered [3, [l^. In the case of contact interac- 
tions, a general renormalization group analysis of elastic 
one channel pot ential scattering has been studied previ- 
ously in Ref. |16j in momentum space within a Lippmann- 
Schwinger framework using a sharp momentum cut-off, 
which separates between the low and high energy region. 
The standard low momentum expansions for both the re- 
action matrix or the K-matrix (in fact the effective range 
expansion) arise as energy perturbations around the triv- 
ial and nontrivial scattering fixed points. The situation 
gets more involved when long distance forces are included 
and has been tackled in Ref. providing the EFT un- 
derstanding of the long distance modified effective range 
expansion 18]. This modified effective range expansion 
has been used to eliminate the one pion exchange (OPE) 
effects in the ^5*0 channel in Ref. The method of 

Ref. flf\ has been applied more recently for the study 
of peripheral waves [20] and is based on delta shell reg- 
ularization in coordinate space for the short range part 
of the interaction, but cutting off the large momentum 
components. This implies, in particular that the long 
range piece extends down to the origin. In the case of a 
singular potential at the origin this procedure becomes 
ill defined, because a sharp cut-off in momentum space 
does not suppress the short distance components entirely, 
unless states of high angular momentum with a short 



distance suppression cancelling the singularity are con- 
sidered 20] . The relation of Wilson like renormalization 
and power counting has been treated in Ref. [2l| in the 
particular case of the singular tensor component to the 
OPE potential (see also Ref. [H [Hill and [H, HI] for 
the two pion exchange (TPE) extension.). A momentum 
space treatment of Wilsonian renormalization ideas has 
also been proposed in Ref. [27i] allowing for a determi- 
nation of model independent low momentum pot entials 
out of several realistic NN potentials [H, [2^ l30l ] . The 
role of redundant operators in the absence of long range 
potentials has been discussed in Ref. [sij . 

A particular advantage of the BC has to do with the 
treatment of non-perturbative renormalization when long 
distance potentials are included. Although the momen- 
tum space treatment accommodates more general situ- 
ations such as non-local potentials than those described 
here, we feel that for the most frequent case of long dis- 
tance physics with local potentials the analysis in coor- 
dinate space becomes more transparent ^. Another im- 
portant reason to prefer coordinate space in our analysis 
is that the Schrodinger equation defines a second order 
boundary value problem, and hence a sharp separation 
between long and short range physics becomes very nat- 
ural. In particular, the short distance unknown physics 
may be handled in the spirit of old and modern works 
as a general boundary condition on the wave function at 
the origin. The question of how the origin should be ap- 
proached is delicate, and depends on the postulated long 
range potential. We will discuss this issue along this pa- 
per in detail. 

The paper is organized as follows. In Sect. jlTI we elab- 
orate on the boundary condition regularization, and its 
advantages as compared to other coordinate and momen- 
tum space regularizations. In Sect. (jIII[) we make a com- 
prehensive renormalization group analysis of boundary 
conditions. The corresponding infrared fixed points are 
determined and identified, and the renormalization group 
fiow due to energy perturbations and potential perturba- 
tions is discussed. For power-like singular potentials at 
the origin we also establish the fixed points, limit cycles 
and attractors. In Sect. (jIV[) we establish the relation be- 
tween the renormalization group fiow of boundary con- 
ditions and the well known variable phase approach to 
potential scattering. This provides a nice interpretation 
of the boundary condition regularization, which suggests 
several working schemes. The generalization of higher 
partial waves and coupled inelastic channels is rather 
straightforward and is presented in a sketchy manner in 
Sect. (|V|. Finally, in Sect. (|VI[) we draw our conclusions. 

Since the present paper had its main motivation in 
the study of the NN interaction problem in effective field 
theories, most of the examples considered along this work 



^ Actually particle exchange implies in the non-relativistic limit a 
local and perhaps energy dependent potential 
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are taken from NN scattering in the ^Sq singlet channel. 
According to Weinberg's power counting, the long range 
reduced potential between the nucleons in this channel 
can be written as a low energy expansion which takes 
the form 0, HI 

Mn Vnn = C^NN = C^LO + C^NLO + C^NNLO + ■ • • (3) 

where LO refers to leading order, NLO to next to lead- 
ing order, NNLO to next to next to leading order, and 
so on; J/lo is the well-known one pion exchange (OPE) 
potential, while C/nlo and J/nnlo come mainly from two 
pion exchange (TPE), although they contain some mi- 
nor contributions to the OPE potential. In this work, for 
convenience, we will mainly use the terms OPE and TPE 
for the potentials, which refer to 

f^OPE = Ulo UtpE — Ulo + f^NLO + f^NNLO • (4) 

It should be noted that in the ^5*0 channel the OPE piece 
is just the well known Yukawa potential 

UoPEiR) = -^e-"^'', (5) 
an 

where a = I&tt f"^ / [M ^ g^) = 0.7 fm and m — 
0.699 fm~^ is the pion mass (/ — 0.468 fm~^ is the pion 
weak decay constant, g — 1.29 is the pion axial coupling 
and Mm = 4.758 fm"^ is the nucleon mass). For the TPE 
potential, it is enough to know that for distances below 
the pion Compton wave length, mR <C 1, it behaves as 

UtMR)^-^. (6) 

where a ~ 1.64 fm in the singlet channel, although it 
depends on the set of parameters one uses (for further 
details on the conventions used along this work for the 
TPE potential, see Refs. [li,!!!]). 

II. BOUNDARY CONDITION 
REGULARIZATION 

For simplicity, let us consider the Schrodinger equa- 
tion for S-wave scattering with reduced potential U (r) = 
2^V{r) and reduced wave function u(r), 

-ul(r) + U{r)uu{r)=euk, (7) 

with the asymptotic behaviour at infinity 

Uk{r) ^ sixv{kr + 5{k)) , (8) 

and subject to the mixed boundary condition 

u'^{R)~-Lk{R)uk{R) = {). (9) 

The coefficient Lk{R) encodes the physics below the scale 
r < R down to the origin at the momentum k ^. For a 



^ This is the most general boundary condition which makes the 
Hamiltonian self-adjoint in the interval R < r < oo. The hard 



given value of the inner boundary radius R, we get a solu- 
tion u{r, R) which depends both on the distance r and on 
the inner boundary radius R. Obviously, the phase shift 
d{k) inherits this R dependence. The boundary condi- 
tion represents our lack of explicit knowledge at some 
low scales, r < R, while we assume complete information 
on the potential, actually a local one, U{r) for r > R. 
For a given R we do not expect the scattering phase shift 
to depend strongly on this lack of explicit knowledge for 
wavelengths larger than R, i.e. kR ^ 1. We will show 
below how this statement can be made more quantita- 
tive. 

Finally, as in any method we want ultimately to remove 
the regulator, i.e. to take the limit of the boundary radius 
to zero, R —^ 0. This implicitly requires to extend the 
potential to short distances. In practice, we expect that 
for a radius much smaller than any other length scale 
in the problem a smooth limit should be obtained. As 
we will see in practical cases of interest in NN scattering 
some cautions must be taken, because although the limit 
in the bare parameters is not necessarily smooth, the 
physical results turn out to be indeed well behaved under 
certain circumstances. 

The BC regularization method has several advantages 
over other methods. For clarity, we list them here, 
although their complete meaning will become obvious 
along the paper. 

• It does not break any symmetry. In the present 
context this is a trivial statement, because it is ap- 
plied to potential scattering once the CM motion 
and angular dependence are separated. 

• It can incorporate any higher order derivative in- 
teractions. This is easily done by making an en- 
ergy expansion of the boundary condition. This 
excludes the subtraction method of Ref. [s^ . 

• It is a non-perturbative regularization. This is par- 
ticularly interesting if one wants to discuss pertur- 
bative approximations, or power-counting schemes. 
A good test to those perturbative treatments is 
to see whether higher order corrections are indeed 
small by comparing with the full non-perturbative 
solution. In dimensional regularization (DR) there 
is no way to make such calculations, unless the reg- 
ulators are analytically removed [SJ) 11^ or the po- 
tential has a particularly simple separable struc- 
ture d^. 

• It can be applied to attractive or repulsive sin- 
gular potentials at the origin A crucial fea- 



core boundary condition, u(R) = corresponds to formally take 
the limit L — > oo needs a separate discussion, since there is no 
continuous dependence on the parameter L at L = oo. 
^ By singular potential we understand a potential fulfilling the 
condition lim, — ,()r^\U{r)\ = oo or lim^^o ^■^'^('') < —1/4. 
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ture for these potentials is that they are non- 
perturbatively renormahzable but become pertur- 
batively non-renormalizable [s^ . Again, dimen- 
sional regularization both in the minimal subtrac- 
tion (MS) or power divergence subtraction (PDS) 
scheme is unable to handle this problem, and to 
date there is no calculation dealing with singular 
potentials in DR. Also the delta shell regularization 
in coordinate space at short distances is excluded, 
because it assumes the singular potential to act at 
distances below the short distance regulator [17j |. 



• It generates a one valued renormalization group 
flow because it involves one distance scale only. In 
contrast, regularization by a potential depends on 
at least two distance scales: the range of the poten- 
tial as well as the strength {1/^/U has dimensions of 
length) of the potential. The square-well short dis- 
tance regulators advocated in Ref. ^3§, '3^ provide 
a multiple branched RG flow structure, which does 
influence the phase-shifts. It is not clear from that 
work which branch should one take, a priori. In 
fact, these infinite branches reflect the well known 
non-uniqueness of the inverse scattering problem, 
rather than the short distance singularity of the 
long range potential. 



• It allows a numerical elimination of the cut-off. In 
the calculations we will show below, we will always 
make sure that physical results are insensitive to 
the actual value of the short distance regulator. 
This improves e.g. on Ref. [40] where a finite cut-off 
was imposed. 



• It provides a one-to-one mapping between the ex- 
pansion of the BC and the physical scattering am- 
plitude or inverse amplitude, both for the case of 
natural and unnatural scattering length. In other 
words, if the short distance wave function is trun- 
cated to a finite given order in a low energy expan- 
sion, the amplitude is also truncated at the same 
order. This property is shared with dimensional 
regularization in the MS scheme in the case of small 
scattering length. 



• It is a uniquely defined regularization in the sense 
that the logarithmic derivative of the wave func- 
tion computed from the asymptotic one can be 
uniquely determined from experiment by just in- 
tegrating Schrodinger equation from infinity down- 
wards to the origin. 



some of the results to come, we propose to solve ^ 

- <(r) + U{r)uk{r) - k\k{r) , (10) 

subject to the boundary condition at the origin and nor- 
malization at infinity 



<(0+) 
Uk{Q+) 

Uk{r) 



= k coi 5 s{k), 

sin(fcr + (5(fc)) 
sin (5(fc) 



(11) 
(12) 



We use the notation Ufc(0+) = limj^^^o+ Uk{Rs), since 
we will see that in general a limit must be taken. In the 
absence of long range potential U{r) =0 the phase shift 
is given by 6s{k). On the other hand, if we take Ss{k) — 
we get a standard problem with a regular boundary 
condition at the origin, Uk(0) = 0. The actual problem 
is that Ss{k) is unknown. At low energies both the full 
phase-shift S{k) and the short distance phase-shift Ss{k) 
can be described by some low energy approximation, like 
e.g., an effective range expansion. 



k cot Ss{k) = 
k cot S{k) — 



1 



1 



(13) 
(14) 



where ao.s is the zero range scattering length, ro^s the 
zero range effective range, and ao and tq the full ones. If 
we also make an expansion at low energies of the reduced 
wave function 

Uk{r)=uo{r) + k^U2{r) + ... (15) 

we get a recurrent hierarchy of equations, namely 



~u'iir) + U{r)uoir) = 0, 
as<(0+) + uo(0+) - 0, 
uo(r) 1 - 



(16) 



r 
a 



at zeroth order and 

-~u'^{r) + U{r)u2{r) 
a5W2(0+)+«2(0+) 
U2{r) 



uoir) , 

^rsasUo{0+) , 



(17) 



6a 



Bar -I- 3aro) , 



at second order and so on. These equations suggest 
a scheme to proceed in practice. If the short distance 
physics could be deduced entirely from the potential we 
would set Q!o,s = 0, ro,s = 0; and so on. Then, the full 
phase shift 5{k) and hence the full low energy threshold 



One of the problems which this regularization manages 
to deal with rather transparently is the disentanglement 
between short and long distance physics. Anticipating 



For the time being we will assume a completely regular potential. 
Singular potentials will be discussed separately below. 
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parameters would be determined entirely from the solu- 
tions of the regular problem at the origin ^. On a leading 
order (LO) approximation, one can improve on that by 
treating a and U{r) as independent variables and predict 
6 and the remaining parameters of the effective range ex- 
pansion , i.e. ro, V2, and so on. In the next-to- leading 
order approximation (NLO) ag, and U{r) are regarded 
as independent variables. 

The standard way to proceed would be to integrate 
the equations, Eq. (|16|) . (|17p and so on, from the origin 
(or a sufhciently small radius Rs) and then to adjust the 
short distance parameters to get the proper threshold pa- 
rameters. Instead, one can simply integrate from infinity 
downwards, with a known value of ag, using Eq. (|16p to 
obtain Q!o,s and then one can use Eq. pO|) together with 
Eq. and Eq. (fT5)) to compute 6{k) for any energy 
with a given truncated boundary condition. This proce- 
dure provides by definition the low energy parameters we 
started with and takes into account that the long range 
potential determines the form of the wave function at 
long distances. The only parameter in the procedure is 
the short distance radius Rs, which we expect to pro- 
duce a smooth limit for the phase shift when we remove 
it by taking the limit Rs (in practice Rs should be 
smaller than any other length scale in the problem). 

Unfortunately, numerical downwards integration of the 
Schrodinger equation is a rather unstable and delicate 
procedure because at short distances the irregular so- 
lution starts dominating and high precision may be re- 
quired to determine the low energy parameters at short 
distances, in a way as to recover the long distance wave 
function. This causes a sort of practical and spurious 
irreversibility triggered by the irregular solution; down- 
wards and upwards integration may not necessarily be 
faithfully represented as inverse operations of each other 
at the numerical level (see e.g. [l^ for further details). 



III. RENORMALIZATION GROUP ANALYSIS 
OF BOUNDARY CONDITIONS 

A. RG equation 

We want to determine the evolution of the boundary 
condition, Lk{R), on the boundary radius, R, by impos- 
ing the physical requirement of independence of phase 
shifts. This is very much in the spirit of the derivation 
of the Callan-Symanzik equation as applied to the renor- 
malization of Green's functions in Quantum Field The- 
ory. The purpose is to establish contact with the methods 
of Ref. [ig] where the analysis is carried out entirely in 
momentum space, and to show that the whole discussion 
can quite naturally be carried out in coordinate space. In 
Appendix|E]we show how the method can also be applied 



to the case where one uses a square well potential 39] to 
regulate the short distance physics yielding a multivalued 
evolution. (See also the discussion in Ref. [41].). 

In order to proceed further, we make the infinitesimal 
change of the boundary radius R R + AR and take 
into account the total derivative 



du{r, R) 
OR 



UR(r, R) 



(18) 



Then, the derivative of the boundary condition with re- 
spect to the boundary radius is given by 

u"{R,R) + u'r{R,K)- L'k{R)u{R,R) (19) 
- Lk{R){u'{R,R) + UR{R,R)) = Q. 

Deriving also Schrodinger's equation with respect to the 
inner boundary radius R we get 

- <(r, R) + U{r)uR{r, R) = ^^^^(r, R) , (20) 

and the asymptotic wave function 

w(r, i?) -> sin(fcr -f (5(A:)) , (21) 

we get 

un{r,R) cos{kr + S{k))SR{k) , 
u'{r,R) cos{kr + d{k)) , (22) 
u'j^{r,R) - sm{kr + S{k)) SnXk) . 

Thus, using Lagrange's identity we get 

— —uru" — u'^u = {—Ufju' + u'j^u)' . (23) 

Integrating between R and oo and using the boundary 
condition, Eq. ^ and Eq. (f23|) . we finally get 

-k^= [k' - U{R) + 4(i?) + Lk{Kf] u{R, Rf . 

(24) 

This equation tells us how the phase shift changes as the 
inner radius is changed. If we require the phase shift not 
to be dependent on the particular choice of R (renormal- 
ization group invariance) we get ^ 



LiiR) 



U{R) + Lk{Rf 



(25) 



This equation governs the evolution of the boundary con- 
dition, i.e., the logarithmic derivative of the wave func- 
tion at the boundary, which shows that in order to guar- 
antee independence of the phase shifts with respect to R 



We remind that we are assuming a regular potential at the origin. 



® This defines RG invariance by a continuous change in the cut-off 
parameter. Note that there is another possible solution to the 
equation, namely the zeros of the wave function, Uj.(_R„(fc), R) = 
0, at some discrete set of short distance cut-offs Rn{k) which ob- 
viously depend on energy. This is equivalent a hard core poten- 
tial, as opposed to the soft boundary condition (|9j. We analyze 
this possibility later on. 
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at all energies the boundary condition must also depend 
on energy, a not surprising result. 

It is instructive to see that Eq. (|25p can be almost 
trivially deduced, by applying Schrodinger's equation to 
the boundary condition at the point i?+ , but the relation 
to renormalization and the approach of Rcf. 16] is less 
obvious. 

Eq. (|24p has interesting consequences as regards the 
low energy limit of the boundary condition and provides 
in addition an error estimate of the phase shift for fc ^ 0. 
If we go to the zero energy limit, the phase shift behaves 
as (5 ~ —ak, with a the s-wave scattering length and 
Eq. ([24|) becomes, 



da 
dR 



U{R) + L'^{R) 
u{R, Rf 



L^{Rf 



X lim 



fc2 



(26) 



The limit in Eq. (|26p is finite since the normalization at 
infinity of the wave function in the limit fc — + is given by 
M ~ fc(r — a). The independence of the scattering length 
with respect to the BC implies the energy independent 
evolution equation 



L'^{R)^-U{R)+L^{Rf 



If we assume this equation we get 



1 d5 
kdR 



k 



2u{R, Rf 



(27) 



(28) 



where we have used that in our normalization u = 0{k) 
for small k. This means that by making the physical 
scattering length BC independent the phase shifts are R 
independent only to order fc'^, if i? 1/fc. This argument 
can be extended to higher orders in kR; if we solve the 
evolution equation to order fc" the error in the phase shift 
is C'(A;2"+2). 

The previous estimate has, in addition a direct appli- 
cation in the renormalization of singular potentials (both 
attractive and repulsive) at the origin, of the form 



U{R) 



1 / a 

Vr 



(R^O). 



(29) 



for n > 2. If we approach the origin R —^ 0, we may 
neglect the energy term in Eq. p4|) . This means that 
the condition for the phase shift to be independent on 
the boundary condition becomes the condition of inde- 
pendence of the scattering length. This, in turn means 
that if the theory is renormalizable at zero energy it is 
renormalizable at any energy. We will see below an al- 
ternative and more appealing formulation of this fact. It 
has been known for a long time by a detailed study of 
the wave functions close to the origin using the WKB 
approximation (37j . and it is very rewarding to provide 
such a simple derivation of this result within the present 
framework. 



The BC parameter Lk{R) has dimensions of inverse 
length, so it is natural to measure it in units of the bound- 
ary radius i?. 



Lk{R) 



UR) 
R 



The equation satisfied by [R) is 

R^^iu{l-ik) + [U{R)-e] R' 



(30) 



(31) 



This is a Ricatti type equation. By using the superposi- 
tion principle for the wave function 

Uk{R) ^UkAR) + k cot S{k)ukAR), (32) 
with Uk,c — > cos(fci?) and Uk^s ^ sm{kR)/k we have that 



ik{R) = R 



'kiR) _ ^<jR) + kcot5ik)u'^JR) 



= R 



Uk{R) Uk,c{R) + k(^otS{k)uk,s{R) 



(33) 



whence the phase shift can be explicitly determined. Us- 
ing the independence of the phase shift on the cut-off ra- 
dius it can be easily shown that two solutions S,k{R) and 
^k{Ro) are related by a Moebius bilinear transformation 
(see Ref. (H), 

^ /n\ A(i?, fio)^fc(fio) + B{R, Rq) 

CiR,Ro)URo) + DiR,Ro)' ^ ' 

where A, B, C and D depend on the potential U only. 
The former expression disentangles the boundary condi- 
tion parameterizing the unknown short distance informa- 
tion from the known long range potential. Actually, the 
matrix 

M(i?, Ro) - (^^(^^ ^(^^ j (35) 

satisfies the group properties 

M{R, R')M{R', R") = M(i?, R") , (36) 

which faithfully represents the dilatation group for the 
short distance cut-off R XR. We will study below 
the stability structure of this group corresponding to the 
the infrared limit i? — > oo and the ultraviolet limit R — > 
0. For a study of the periodic case in log(i?) see e.g. 
Ref. fH. 



B. Long distance Fixed Points at zero energy 

For a potential of typical size a and at low energies 
ka <^ 1 we may look at the region a <C i? <C 1/fc where 
both the potential and the energy can be neglected, yield- 
ing the equation 



a<i?<l//c. (37) 
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Note that the equation is scale invariant under the change 
R XR within the interval a <^ R <^ 1/k. Obviously, 
having a finite interval breaks the scale invariance, since 
the interval boundaries also change. Thus, both the po- 
tential and the energy are scaling violating perturbations. 

Fixed points in the dimensionless boundary condition ^ 
are defined as those fulfilling the condition = 0. Hence, 
Eq. ([57)1 has two fixed points at = and = 1- The 
point ^0 = is unstable since any small perturbation 
of its value at say r ~ a results in increasingly large 
deviations from the fixed point, as can be seen directly 
by analyzing the differential equation, Eq. ([57)1 . On the 
contrary, the point = 1 is stable. Since both fixed 
points are associated with long distance behaviour they 
correspond to infrared (IR) fixed points. The physical 
interpretation of fixed points in the present context is 
clear; by varying the BC in the range a ^ R ^ 1/k 
according to Eq. (|37p , we guarantee physics independence 
at low energy. The fixed points represent very special 
situations where the BC itself does not change with the 
boundary radius, i.e. it is scale invariant. For such a 
situation, we can characterize the boundary condition by 
a number instead of a function of the boundary radius. 
In a real case we do not expect to have exact fixed points, 
but only some approximation to them. For this case it 
will be interesting to choose the scale a ^ i? <C 1/fc 
to describe the boundary condition function because we 
expect a weak dependence on the scale. 

For the unstable point, = Oj the situation is such 
that an extreme fine tuning of the BC is required to have 
this scale invariance. On the contrary, the stable fixed 
point, 5o = 1; does not require this accurate fine tuning. 

To analyze the physical situation corresponding to 
these fixed points let us now take into account that 
Lq(R) = u'q{R)/uo{R). In the region a < i? < 1/fc 
we have the asymptotic wave function (we use the nor- 
malization condition ito(O) = 1), 



which solution is 



uoiR) 



1 



R 
a 



yielding 



Lo{R) 



R 



R-a 



(38) 



(39) 



Thus, the unstable fixed point ■Co = corresponds to a ^ 
00 (zero energy bound state) and the stable fixed point 
^ = 1 to a = (trivial scattering). These conclusions are 
in full agreement with the momentum space analysis of 
Ref. M- 



C. Positive Energy Perturbation 

In the region where the potential does not act r 3> a, 
the equation satisfied by ^k {R) is 



R^=Ul-^k)-k'R'.. 



(40) 



^kiR) ^ kRcot [kR + 6{k)] 



(41) 



We have fixed the arbitrary constant by imposing that 
the extrapolation of the logarithmic derivative to the ori- 
gin can be related to the scattering phase shift, S{k). 
Then we have 



k cot d{k) = k 



kR + ^k{R) cot{kR) 
kRcot{kR) - ^kiR) 



(42) 



Thus, an expansion of the scaled boundary condition in 
the form 



(43) 



valid in the region a <C i? <C 1/fc yields an expansion for 
the phase shift as given by Eq. (|^T|l . 



k cot 6{k) 



1 & , ^o'- 3^0 + 36 



Rl-^o 



m - 1)^ 



(44) 



Whereas written in this way the expansion for k cot S 
around the trivial fixed point ^ = 1 induces increas- 
ingly large contributions for increasing orders in k, spoil- 
ing the convergence of the low energy expansion, in the 
case of the nontrivial fixed point ^ = this yields a 
perfectly well defined expansion. A similar situation oc- 
curs for the expansion in tanS/k, although with opposite 
fixed points corresponding to the divergent and conver- 
gent case. Thus the perturbation theory around the non- 
trivial fixed point (large scattering length) corresponds to 
an effective range expansion of the form 



k cot S{k) 



(45) 



whereas for the trivial fixed point (small scattering 
length) the low energy expansion reads. 



tan(5(fc) 



-a - (3k^ - jk'^ + 



(46) 



Matching both expansions. Eq. and Eq. we get 
the identifications 



« 1 2 



(47) 
(48) 



The convergence radius of these low energy expansions 
has to do with the longest distance singularities of the 
potential. For Yukawa like behaviour, U ^ e^^^/i?, one 
has |fc| < to/2, due to the branch cut at fc = ±zto/2. 
Comparing the expansions (|43p and we have 
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Co(i?) 



R 

R-a ' 

QaR^ - 2R^ + (-6i? + 3ro) 

6 (a - Rf R ' 
- (-24ai?5 _^4i?6 - 15a3 (4i?3 - 6 i?2ro + 3 i?rg - 12^2) + iQa^ {iR'^-R^r^ - 6 RV2)) 

180 R^ {-a + Rf 

I 



(49) 
(50) 

(51) 



w 






:'; \ \ 1 ■ 
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--'>' ' / ii 










: / / 
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-0.5 




^0 



0.5 



FIG. 1: The relation between the s-wave parameters ^2 and ^0 
for the pure short range theory defined as the energy dimen- 
sionless coeflicients of the logarithmic derivative of the asymp- 
totic wave function Ru'^iR) /uk{R) = £,o[R) + {kRfS_2{R) + 
. . . depending on particular choices of the ratio ro/a between 
the effective range and the scattering length. We also plot the 
flow for NN-scattering the case of singlet ^So and triplet '^Si 
states. 



If we express the energy correction to the scaled bound- 
ary condition ^2 in terms of ^0 by eliminating the distance 
R we get 



- 3^0 + 3)(3ro/a - 2) - 3ro/a 
6^0 



(52) 



The renormalization group flow corresponding to the de- 
pendence of ^2 in terms of ^0 is depicted in Fig. ((T|) for 
several values of the ratio between the effective range rg 
and the scattering length a. Similar pictures could be 
obtained for any other pair of ^„ variables. For illus- 
tration purposes we also show the flow for the case of 
singlet ^5*0 and triplet ^S'l states in the NN interaction 
where ao = — 23.74 fm, tq = 2.77 fm and ao — 5.42 fm , 
rg = 1.75 fm respectively. 



D. Zero energy Short distance critical points 

1. Fixed points and Cycles. 

We analyze now the short distance behaviour of 
Ea. (j3ip corresponding to the ultraviolet regime. Gen- 
erally speaking, if the solution to the Schrodinger equa- 
tion is written as a linear combination of a regular and 
irregular solution, uo{r) — cq Uirrcg('') + Ci Mrcg(f ), then 



Co(i?) = R 



(i?) -I- Cl-Urcg(i?) ' 



(i?) + cim' (i?) 

H- 



(53) 



where the ratio cq/ci can be fixed by choosing ^o(^o) 
with Rq some reference scale. Obviously, regular solu- 
tions are unstable fixed points whereas irregular solutions 
correspond to stable fixed points since Uincg always takes 
over when R —^ 0. 

We will first consider (Sect. IIIID2P the case of regular 
potentials, for which the following condition is fulfilled 



lim i?^ U(R) = 



(54) 



and afterwards the case of singular potentials which be- 
have as 1/i?" at short distances. In this latter case we are 
going to make a further distinction between those power- 
law potentials with n > 2 (Sect. IIII D 3[) and the inverse 
square potential, U{R) = g/R^ (Sect. IIIID4p . For a 
short review of the short-distance solutions of singular 
power-law potentials see Appendix \^ 



2. Regular potentials. 

For a regular potential the wave function behaves lin- 
early at short distances, u{R) = cq + ci i?, yielding 



ciR 



Co 



(55) 



which has the fixed points ^0 = 0, 1. The first one, 
^0 = is stable and corresponds to irregular solution 
of the Schrodinger equation. The convergence towards 
this fixed point is linear in R 



UR) 



Cl 

Co 



R + 0{R^ 



(56) 
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The other fixed point, = 1 is unstable and corresponds 
to the regular solution u{R) ~ R. 

A special case is that of the Yukawa potential 



a R 



(57) 



which has al/R singularity for short distances. For this 
potential the behaviour of the wave function at short dis- 
tances is given by 



u{R) = Co 



1 + mR 



3R i?, ,R, 

log(-) 

2a a a 



+ ciR, (58) 



for which the same fixed points than discussed before 
for a common regular potential are reproduced (namely 
^0 = 0, 1) , but with logarithmic convergence towards the 
stable fixed point = 



Co(i?) 



Cl 

Co 



R + mR- 



^-^-^log{^)+0{R'),{59) 
la a a 



from which an R logi? convergence trend is deduced 
(which, incidentally is independent of the actual ci/co 
ratio). 



choosing ^o(-Ro)- Close to the origin we can use the 
asymptotic expansions 



TTzcos ((i//2 + l/4)7r — z) (63) 

and hence we have 

eo(i?o)cot$-l 



eo(i?o)+cot$ ' 



with 



n - 2 



R 

Rq 



l-n/2 



(64) 



(65) 



where Rq is some short distance reference scale. So we 
see that for n > 2 we have an attractor {R,^o{R)) with 
a fractal asymptotic {R — > 0) dimension d = 2 — 2/n. 
The fractal dimension can be deduced from the scaling 
properties of the zeros of ^q{R), R^'^~'^ ~ 1/-^, being 
Rn the N-th zero. 



4- Inverse square potential 



3. Power-law singular potentials. 



For a power-law potential which behaves as 



U{R) 



1 / a\" 
\R 



n>2. 



(60) 



as i? — > we have one scale only, so we can define the vari- 
able i? = ax, in such a way that the coupling constant 
becomes one. The solution to the Schrodinger equation 
for the repulsive case reads 



u{x) = Cq^/xK_ 



X 



l-n/2 



— \n/2- 1 

2,l-n/2 



^ \nl2 - 1 J ' 



(61) 



where co,i are integration constants and Ki^{z) and I,j{z) 
are regular (exponentially increasing) and irregular (ex- 
ponentially decreasing) modified Bessel functions respec- 
tively. The fixed points are ^o = ±oo corresponding to 
take cq — Q and ci = respectively; ^o = +oo is unstable, 
while ^0 = ~oo is stable. 

For the attractive n > 2 case, the solution to the 
Schrodinger equation for the attractive case reads 



u{x) 



C\\fxJ_ 



X 



l-n/2 



^ \n/2- 1 

^l-n/2 



+ 



n/2~l 



(62) 



where are oscillating spherical Bessel functions. From 
here one can obtain the dimensionless logarithmic deriva- 
tive S.oiR), while the ratio cq/ci may be determined by 



The inverse square potential, U — g/R^, requires a 
separate study. Note that in that case the RG equa- 
tion, Ea. (j3ip . is formally invariant under the continu- 
ous scaling transformation, R —^ fiR. This symmetry is, 
however, explicitly broken by an initial condition ^o(-Ro)- 
Using Eq. (pij) we get the zero energy fixed points at 



(66) 



for g > —1/4. The positive and negative roots corre- 
spond to stable and unstable fixed points respectively. 
For g < —1/4 we have purely imaginary solutions. This 
is the signal for a limit cycle. Actually, the solutions are 
given by 

uo{R) = coR^- +ciR^+ 



A± = (l±yr+4^)/2 for 5>-l/4, 

(67) 



(68) 



ua{R) = cqV Rcos{XlogR) + ciV Rsm{X\ogR) 
A = a/-1 - ig for g < - 1/4 . 

So, in the first case we obtain 



^oiR) 



coX-R 



A- 



ciX+R^+ 



CqR^- + Cli?^+ 



for g > -1/4, (69) 



which has an attractive fixed point at ^o = ^-i a-nd a 
repulsive one at ^o = A+. In the latter case, g < —1/4, 
we obtain 



Co(i?) = A cot 



-1 2A , , ^ 

2UR^l+''^^Ro^ 



1 

2 ' 

(70) 
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FIG. 2: Running of ^o(^) for regular and singular potentials. In the left panel we show ^o(^) for the case with no potential in 
which £,oiR) = R/(R — ao), for a square well potential and for a Yukawa potential (the One Pion Exchange potential in the ^So 
channel). In the right channel we show (,o{R) for the case of the Chiral Two Pion Exchange potential which for short distances 
behaves as —a^/R^; the lines representing the renormalization flow of S,o{R) becomes more dense in its way to the origin, giving 
rise to a nontrivial fractal dimension d = 2 — 2/n = 5/3in the i? — > limit. These examples are taken from neutron-proton 
scattering in the singlet s-wave channel, in which the scattering length is qq ~ — 23.74 fm, and R is expressed in fm. 




■4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 -4 -2 2 4 6 8 10 0.1 0.2 0.3 0.4 0.5 



FIG. 3: Running of vs ^2 for regular and singular potentials (note the log scale). In the left panel we show the case of a 
Yukawa potential (the One Pion Exchange potential in the ^So channel). In the middle channel we show the case of the Chiral 
Two Pion Exchange potential. In the left panel we compare the case with no potential to that of a OPE and TPE potential. 
These examples are taken from neutron-proton scattering in the singlet s-wave ^So channel, in which the scattering length is 
ao — — 23.74 fm, and ro = 2.77fm. These plots should be compared with Fig. ([l}. 



In this case the discrete scahng property ^o{Re^'"^^) = 
^o{R) with N = 0,±1,±2, ... typical of the hmit cy- 
cles ^1] is reproduced. For instance, if we consider 
£,o{Rn) — 0, then the sequence of the zeros of is given 
by Rn+1 = e^'/^RN- This corresponds to the Russian 
doll renormalization (see e.g. Ref. [4^). Coordinate and 
momentum space analyses have been treated in [46] and 
[irf respectively. 



5. Overview of the ultraviolet limit 

On the light of the previous discussion we obtain 
for the short distance dimensionless boundary condition. 



^o(^^), the following behaviour 

• For a regular potential we have two fixed points, an 
attractive one, corresponding to the irregular solu- 
tion of the Schrodinger equation, and a repulsive 
one, corresponding to the regular solution. This 
means that all solutions go to the irregular solu- 
tion at the origin. 

• For a singular potential with n = 2 and coupling 
g > —1/4 we have two fixed points. For g < —1/4 
there are limit cycles. 

• For a repulsive singular potential with n > 2 we 
have two fixed points. The attractive one corre- 
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spends to the irregular solution. 

• For an attractive singular potential with ?i > 2 we 
have an attractor with asymptotic fractal dimen- 
sion d = 2 — 2/rt. 

In Fig. ([2]) we present three of the possible situations. 
The examples are taken from neutron-proton scattering 
in the singlet s-wave ^5*0 channel, for which the scat- 
tering length is ao = — 23.74 fm. In the regular poten- 
tials example we use for the discussion the One Pion Ex- 
change (OPE) potential, which, for this channel, takes 
the form of an usual Yukawa potential as the one given 
in Eq. dST]), for which m = 0.699 fm"^(= 138 MeV) 
and a — 0.7 fm, and a Square Well potential with range 
i?o = 1/m ^ 1.43 fm and depth Uo = -0.1 fm^ In the 
singular potential example we use the chiral Two Pion 
Exchange (TPE) potential, which behaves as —a^/R^ at 
distances below the pion Compton wave length. For com- 
pleteness we draw also in Fig. ([3]) the behaviour of vs. 
£.2 very much in spirit of our study of Sect. IlllCl on pos- 
itive energy perturbations (see Fig. ([T]) ). 



E. Low energy expansion of the BC with a 
potential 

If we relax the condition i? 3> a, we can undertake a 
low energy expansion to get the set of differential equa- 
tions for the low energy scaled BC coefficients, 



R 



da 
dR 

^dR 



eo(l - eo) + UiR)R^ , 
-6(1 + 2^0)-!, 



-^4(3 + 2^0)-^: 



2 ; 



(71) 
(72) 
(73) 



valid for any R. The r.h.s. are the corresponding beta 
functions of the renormalization group flow. Note that 
unlike the standard RG equations a manifest scale de- 
pendence shows up due to the potential U{R) which de- 
scribes the long range physics. The fulfillment of these 
equations guarantees cut-off independence of low energy 
parameters, a, rg, V2 and so on for any value of the cut- 
off Rs- Moreover, these equations exhibit a natural hi- 
erarchy; the solution of a given coefficient ^2n depends 
only on the previous ones £,2n-2, • • ■ ^- Furthermore, 



they are non-perturbative in the potential and do not 
require any off-shell information, as in the momentum 
space treatments • Actually, the solution for the first 
equation in the absence of a potential reads, 

UR) = ^^, (74) 
R — a 

where a is the integration constant which can be identi- 
fied with the physical scattering length. In the presence 
of the potential U (R) , this suggests a solution of the form 



R 



R-aoiR) 



(75) 



where a{R) is an undetermined coefficient, satisfying 
Eq. (|122p of Sect. |TVl We will show in the next section 
that ao{R) is the scattering length corresponding to the 
truncated potential U{r) for r < i? and zero otherwise. 
Actually, making use of the analogy we can solve the 
equations, fTT]) . ((72| and (|73p in a more efficient manner, 
see Eqs. (fm)) . (fT25l) and (fTM]) . 

The set of equations (HIl), ^ and ^ have to be 
solved with some initial conditions. For asymptotically 
large distances we must have ^'q ^ , and hence 



Co(oo) 
?o(oo) 



1 

= 



(a 7^ oo) , 
(a = oo) . 



(76) 
(77) 



Thus, unless a = cxd, all solutions go asymptotically to 
the stable fixed point, according to the general theory. If 
a 3> a then there is a region for r ^ a where ^ 1 and 
^0 remains almost constant. Likewise for a <C a we have 
^0^1 for r > a. Finally, if a ~ a we have ^ 1 only 
in the region r a. 

If we solve around the fixed point = 1 we get 



^o(i?) = 1 + i? / U{r)dr. 



(78) 



Going to R 



oo since already at R • 



a 



a we may have 
(79) 



Thus, if a ^ a we have ^o(a) ~ 1 and if a ^ a then 
^o(a) ~ 0, even if £,o{oo) = 1 (unless a — oo, in which 
case we do have ^o(oo) — 0). Thus, taking the scale 
R ^ a seems like a good choice to be as close as possible 
to the fixed point situation ^' 0. 



F. Error estimates 



1. General Considerations 



This is in contrast to momentum space parameterizations of the 
short distance potential via momentum dependent perturbations 
Va{k' ,k) = Co + C2{k'^ + k''^) + ■ ■ ■ where operator mixing occurs. 
In practice this means that if at zeroth order one keeps only the 
counterterm Co its running is fixed by the regularization method 
and a renormalization condition, most naturally fixing the scat- 
tering length a. When the new term C2 is considered, the run- 
ning of the lower order counterterm Co is modified and both 
Co and C2 are intertwined, their values at a given cut-ofi' being 



Exact renormalization group invariance requires the 
knowledge of the complete phase shifts at all energies for 



fixed by a and the effective range. A way out is to consider en- 
ergy dependent perturbations instead of momentum dependent 
ones [la . llTll and the hierarchy that we find is recovered. 
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the fulfillment of Eq. p5l) . but practical computations 
demand the use of a limited amount of physical infor- 
mation as input in order to have predictive power. One 
example is a theory in which we know the value of the 
scattering length, ao, and consequently we want to fix its 
renormalization group flow, which can be obtained from 
the low energy limit of Eq. ([M)) , yielding 



da 
dR 



= lim 



UkiR, R) 



L'oiR) + LoiRf - U{R) 



0. 



(80) 



Thus, ao(R) — aa{Ra) implies the fulfillment of Eq. ([7T|) . 
but not of Eq. (|72p or higher order ones. This theory can 
be shown equivalent to truncating the boundary condi- 
tion (see Appendix |B| in the R limit 



lim Lk(R) 



lim Lo(R) 



with Lq{R) fulfilling the RG equation 

L'oiR) + LoiRf = UiR) 



(81) 



(82) 



subjected to the asymptotic boundary condition 
Lo{R) ^ l/(i? - ao) at i? ^ oo. 

We can try to improve the description of the phase 
shifts by making RG independent not only the scattering 
length ofQ, but also the effective range rg, or even higher 
order parameters of the effective range expansion. In 
case we fix ao, To, W2, • ■ • , Wn, we can write the boundary 
condition as 

Lfe(i?) = Lo + fc2L2 + --- + fc2»L2„ + ... (83) 

where Lq^ L2, . . . , ^271 obey their respective renormal- 
ization group equations ^, while the higher order terms, 
represented by the dots, do not. 



2. Zeroth Order Truncation 

If we truncate the BC to zeroth order in the energy, 
i.e. Eq. is fulfilled while Eq. ([72]) and higher order 



The relation between the truncation of the boundary condition 
and the fixation of low energy scattering observables, i.e. the 
different terms of the effective range expansion, can be made 
clear by rewriting Eq. 12411 as 



dR 



■fccot<5 = [L'^.{R) + lI{R) + P - 
xu(R, Rf , 

which is obtained in just the same way as Eq. J24I I. but changing 
the asymptotic normalization of M(r, R) to 

u{r, R) cos kr + k cot 5{k, R) , 

k 

from this it is obvious that fixing is equivalent to the ful- 
fillment of Eq. {TTJ, fixing ao and ro to fulfilling Eq. l|7H l and 
Eq. II72I 1. and so on. 



ones are not, then the scattering length is independent 
on the short distance cut-off R. The truncation is made 
in order to fix the scattering length of the system, while 
the effective range, the shape parameter and higher order 
terms of the expansion of k cot 6 are determined by both 
the scattering length and the long range potential. 

This situation is equivalent to take Lk{R) = Lo{R), 
from which the cut-off dependence of the phase shift can 
be deduced by making the substitution Lk{R) — > Lo{R) 
in Eq. ([M]), yielding 



dR 



fukiR,R) 



(84) 



For a regular potential, for which we have the behaviour 
Uk (i?, R)/k ^ Co + ci R for small cut-offs, the phase shift 
shows a linear dependence in the cut-off radius which 
slope increases with the momentum fc, so for small cut- 
off radii we have 



dS 
dR 



(85) 



and consequently we find a linear dependence of the 
phase shift with respect to the cut-off radius, A6 ^ R. 

For an attractive singular potential, which behaves as 
U{R) ^ — 1/i?", the situation changes, since the wave 
function is very much suppressed at short distances, 
Uk{R: R)/k ^ R"-/^ as compared to the regular potential 
case. Then we obtain, for small cut-offs, the behaviour 



^^-k^R-/\ 
dR 



(86) 



from which one can deduce a cut-off error A (5 ~ 
and we recover renormalizability, i.e. cut-off indepen- 
dence, in the small cut-off limit 



,. d5 
hm — — 

R^o dR 



0. 



(87) 



This cut-off dependence is depicted in Fig. (jH), in 
which we show the nucleon-nucleon phase shift in the ^So 
channel as a function of the cut-off R for various center- 
of-mass momenta. The phase shift 5 (k, R) is computed 
using an energy independent boundary condition which 
turns out to be equivalent to fixing the scattering length 
to some given value, as for example the experimental one 
ao = — 23.74 fm, for both the one and two pion exchange 
potentials. The OPE potential is just an usual Yukawa 
potential, while the TPE potential behaves as —a^/R^ for 
distances smaller than the pion Compton wavelength. As 
we can see in the figures, in the OPE case the phase shift 
show a linear dependence in R near the origin, while in 
the TPE case the phase shift becomes insensitive to the 
cut-off much earlier. 

Finally, for a repulsive singular potential behaving as 
U{R) ^ 1/^", the phase shifts shows an exponential 
behaviour as the cut-off is removed. This can be under- 
stood from the behaviour of the wave function at short 
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FIG. 4: Cut-ofF dependence of the phase shifts for NN scattering in the ^Sq channel with the OPE and TPE potentials, which 
behave as 1/R and respectively. For computing the phase shift we use a energy independent boundary condition fulfilling 

the equation L'q{R) + Lq{R) — U{R) — 0, which in turns mean that the scattering length is independent on the cut-off (we take 
ao = — 23.74fm). We show the phase shifts for center of mass momenta of kcm = 0.25, 0.5, l.Ofm"^. 



distances (see Appendix [X)) 



UkiR, R) 



^n/2g+(a/i?,)"/=^-V(«/2-l) 



where we have taken the irregular solution since for an 
arbitrary scattering length this is the UV fixed point and 
the regular solution may be discarded. Then, for small 
cut-off 

^ ^ _fc3 g+2(a/i?)"/^- V(«/2-l) ^ (39) 

dR 

from which one can see that the phase shift develops 
a very strong cut-off dependence when one tries to fix 
the scattering length to an arbitrary value, unlike the 
previous cases. 



3. Relation to Orthogonality Constraints 

The truncation of the boundary condition at zeroth or- 
der is equivalent to the fulfillment of orthogonality rela- 
tions between different energy solutions, which was used 
in Refs. [12, H^, HI] to renormalize and obtain model- 
independent predictions for the ^5*0 singlet and '^Si — ^Di 
triplet (deuteron) channels. Orthogonality constraints 
imply that different energy wave functions fulfill the in- 
tegral relation 



Uk{r)uk' (r) dr = 5 {k ~ k') 



(90) 



If we use the corresponding Lagrange identity between 
Mfe and u'j,, we can rewrite the previous orthogonality 
constraint as 

/>oo 

(fc^-fc'^) / Uk{r)uk'{r)dr ^ u kUk' -UkU k' , (91) 



so for k' 7^ /c, and including a short distance cut-off i?, 
we obtain 



u'kiR) ^ u'k'jR) 
Uk{R) Uk'{R) ■ 

Finally, taking fc' = as a reference state 

u'k{R) ^ u'ojR) 
UkiR) uoiR) 



(92) 



(93) 



or, equivalently, Lk{R) = Lq{R), i.e., imposing orthogo- 
nality between different energy solutions is equivalent to 
truncate the boundary condition at zeroth order. 



4- Higher Order Truncations and Inconsistencies 

Naively one can think that a truncated energy expan- 
sion of the boundary condition would lead to a systemat- 
ically more accurate description of the data. In this sec- 
tion we will analyze this problem and show cases where 
this naive expectation is not fulfilled. 

To begin with we can try to truncate the boundary 
condition, for example, to second order in the energy 



Lk{R) ^ La{K) + k" L2{R), 



(94) 



so Eqs. dZH) and ^ are fulfilled while Eq. ^ and 
higher order ones are not. In such a case both the scat- 
tering length and the effective range are cut-off indepen- 
dent, while the phase shift is not. In a similar fashion 
to the previous case, the cut-off dependence of the phase 
shift can be estimated, yielding 



dR 



L2{Rf 



UkiR, R) 



(95) 
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which can be easily computed for regular potentials. 
For such potentials we can take into account that (see 
Sect. El) 



lim Lk{R) 



k cot Ss{k) 
1 1 



as 



■rs 



(96) 



being as and rs the short-range scattering length and 
effective range. Inserting the previous expression on 
Eq. ([95]) . we arrive to the following result 



d6 



lim —:^^k'>{ — 

R^o dR V 2 



rs\^ /ufe(0,0) 



^0. 



(97) 



This means that by adding more terms in the energy 
expansion of the boundary condition, i.e. fixing higher 
order parameters in the effective range expansion, we will 
find the same linear dependence as in the zcroth order 
case of the phase shift on the cut-off R in the ultraviolet 
limit. With each new order in the energy expansion, 
the slope of d{k,R) in the R ^ limit is progressively 
suppressed in energy with respect to the zeroth order 
case, e.g. by a factor of {k rg)^ at second order, k'^ rs vs 
at fourth order, etc. 

One example of this energy suppression is given by the 
OPE potential in the channel, for which rs = 4.46 fm 
(see Sect. lIVBp . In this case, the cut-off dependence is 
in principle smoothened for k < 0.45 fm~^ (^ 88MeV), 
indeed a bit more due to the higher order terms k'^ , k^, ... 
which have been ignored here. The cut-off dependence 
for different energies is depicted in Fig. ([5]), both for the 
case where a and both a and ro parameters are fixed. As 
we can see, the effective cut-off dependence is still a bit 
smoother for two parameters at fc = 0.5 fm~^, and at fc = 
1.0 fm~^ is only slightly worse than with one parameter. 

In the case of attractive singular potentials the renor- 
malization group analysis becomes more complex. The 
first step is the evaluation of L2, which is computed in 
Appendix [Bl yielding for i? ^ 



L2iR) 



uoiR,R) 



■Aro 



(98) 



where A rp is the short distance contribution to the effec- 
tive range ^, and the zero energy wave function is defined 



by 



uq(R, R) = lim 

k^O 



UkiR, R) 



(99) 



The similarity with the L2{R) rs/2 result for regu- 
lar potentials is striking. In fact the formula above can 
be used to obtain rs for a regular potential (for which 
uo{R, R) goes to a constant value at short distances). 

Now we can evaluate the convergence of the phase shift 
with respect to the cut-off for = Lq + fc^ L2, which is 
given by 



dS 
lim 

R^o dR 



Aro 
2 



uq{R, R) 



^0. 



(100) 



This behaviour is depicted in Fig. ^ and at first sight 
looks weird. With every zero of the wave function there 
is a TT jump in the phase shifts, meaning that the sys- 
tem is explicitly sensitive to the appearance of deeply 
bound states, something that should not have practical 
consequences within the domain of applicability of a cor- 
rectly formulated effective theory. There is no obvious 
and unique i? — > limit. 

The previous result seems to suggest that one can- 
not fix both ao and rg and obtain renormalized phase 
shifts. However, a further analysis below will show that 
it is indeed possible and the results are well defined and 
unique. 

Instead of a naive truncation of the boundary condition 
as we have done, we suggest to reorder the expansion as 
follows 



LkiR) = £o(i?) + fc' L2iR) + k^ AiiR) 



(101) 



where Lq and L2 fulfill exactly the RG equations, but 
higher order terms do not (and hence the notation A4, 
Ag, etc, for these terms). In such a case we can see that 
the scale dependence of the phase shift is given by 



^ = -fcS [L2iR)' + 2LoAi + A'^{R)] 



Uk{R, R) 



(102) 



Thus, the total effective range of the system is 

ro = Aro + ro(ao) 

where ro(ao) is the effective range due to the zero energy bound- 
ary condition used to fix the scattering length of the system. It 
can be computed with the well known integral formula 

roc 

ro("o) = 2 / dr [vgir) - ul{r)] , 
JO 

where t)o(r) = 1 — r/ao and i*o(r) is the solution of the zero 
energy reduced Schrodinger equation with the asymptotic nor- 
malization uo{r) vo{r) for large distances. 



In fact, we can consider for example the shape parameter V2, 
which has the following integral definition 

roo 

V2 = dr [vQ V2 - Mo "2] , 
Jo 

where vo(,r) = 1 — r/ao and V2{r) = (r^ — 3 agr + 3 oiq rg) r/6ao , 
while uo and U2 are solutions of Eas. l|16| l and l|17| l respectively, 
subjected to the asymptotic boundary conditions uq — > ^'0 and 
U2 ^ V2- Since the solutions for attractive singular potentials 
behave at short distances as r"/* times some trigonometric func- 
tion, V2 is convergent. A similar argument can be applied to 113, 
V4, and so on, and even to the phase shifts, thus giving a well 
defined limit when the cut-off is removed. 
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FIG. 5: Cut-off dependence of the phase shifts for NN scattering in the ^5*0 channel with the OPE potential when we fix one 
(IC curve) or two (2C curve) parameters, namely the scattering length and the effective range (we take qq = — 23.74 fm and 
ro = 2.77fm). We show the convergence of the phase shifts for center of mass momenta of kcm ~ 0.25, 0.5, 1.0 fm^^ (left, 
central and right panel respectively). We can see that the 2p curve is more convergent than the Ip curve for 0.25 and 0.5 fm~^, 
while only slightly worse for 1.0 fm~^. 




FIG. 6: Cut-off dependence of the phase shifts for NN scattering in the ^So channel with TPE potential (an attractive singular 
potential) when fixing the scattering length and effective range of the system. In the left panel, the boundary condition has 
been truncated at order in the momentum expansion, Lk = Lo + L2, thus giving rise to cycles due to the incorrect 
renormalization procedure used. In the right panel, the boundary condition is expressed as a Fade approximant, such that we 
recover renormalizability in the J? — > limit (we take qq = — 23.74 fm and ro = 2.77 fm). For the convergent case, we show the 
phase shifts for center of mass momenta of kcm = 0.25, 0.5, 1.0 fm^^. 



Then, analyzing the behaviour of Lq, L2 and A4 (see 
Appendix [B]) , we find that 

d S 

— = -k^uo[R,R)u2{R,R) + 0{k^) 
dR 

(103) 

which means that the convergence pattern that emerges 
when we fix both a and ro is similar to that found when 
we fixed ao only, i.e. there is no improvement on the 
short distance scaling suppression. The same arguments, 
but considering Ag, Ag, etc, can be applied to higher 
orders in the momentum expansion of the phase shifts 
yielding identical short distance scaling suppression. 



The nontrivial fact is that fixing both the scattering 
length and the effective range is not exactly equivalent to 
any truncation of the boundary condition Lk{R). All the 
terms in the energy expansion are equally singular, and 
therefore equally relevant. The correct parameterization 
of the truncated boundary condition corresponds to a 
unique Pade approximant looking representation of the 
boundary condition 

L,iR)=<±pk + k^A{R), (104) 

Uq + U2 

where A(i?) is a remainder. Indeed, for short dis- 
tances, the Pade behaves as 1/i?"/^, while the remain- 
der A(i?) ~ R, and hence can be safely ignored when 
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the cut-off is removed. This expression for Lk{R) gives 
an UV scaling for the phase shift, regardless on 

how many low energy parameters a, tq, V2 and so on are 
fixed. The form of this Pade, which could be visualized 
as a renormalization group improvement is driven by the 
Moebius bilinear transformation, Eq. ([34|) which embod- 
ies the dilatation group properties of the short distance 
cut-off given in Eq. p6p . A different Pade approximant 
will lead to spurious cut-off dependences, which will jeop- 
ardize the i? — *■ scaling behaviour. 



5. Remarks on Regular Solutions 

The RG behaviour of regular solutions represents a 
very special case. In it we fix the scattering observables 
to the values corresponding to the regular solution at the 
origin. By doing this the convergence with respect to the 
cut-off is improved noticeably. 

The first case we are going to consider is the triv- 
ial one, that of a regular potential in which we do not 
fix anything, but rather enforce the regular solution 
R) = CqR at a small cut-off radius R. For such 
a case we have the cut-off dependence 



dS 
dR 



k [U{R) 



UkjR, R) 
k 



(105) 



since Lk{R) — l/R. In it i? can be interpreted as the 
starting integration point in any usual integration pro- 
cedure for the Schrodinger equation. For a potential 
which goes to a constant value at the origin, U{R) C/q, 
we have AS ~ R"^, while for a Coulomb- like potential 
at short distances, like the Yukawa potential, we have 
AS ~ i?2 

When one fixes the scattering length to the regular 
value, the convergence rate is given by Eq. ([84]) . In such 
a case, there would be a small improvement for a Yukawa 
potential, which will now converge as A6 ~ R^. These 
results should be compared with the convergence when 
fixing the scattering length to an arbitrary value, which is 
AS R. When one additionally fixes the effective range, 
Eq. ([55)1 describes the convergence rate. In case the ef- 
fective range is taken to be the regular value, we have 
that the short distance effective range is zero, rs = 0; 
L2{R) can be estimated, for example, by the Green func- 
tion methods used in Appendix [Xl giving L2{R) ~ R. 
Then we have AS ~ i?^, in accordance with naive expec- 
tations that increasing the number of counterterms will 
smoothen the cut-off dependence. Further suppressions 
will take place at higher orders. 

The second case is when the potential is repulsive sin- 
gular. As we saw previously, setting the scattering length 
generates an exponentially divergent cut-off dependence 
in the ultraviolet limit, given by Eq. ([55]) . The only way 
in which one can obtain renormalization invariance with 
a repulsive singular potential is by taking the scattering 
length corresponding to the regular solution at the ori- 
gin, for which we obtain a rapidly convergent result for 



the phase shift 

^ ^ i?"/2 g-2(a/fl)"/- V(«/2-l) ^ (106) 

dR 

in which the sign of the exponential has changed with re- 
spect to Eq. ([55)1 , since we are taking the regular solution 
of a 1/i?" potential. When we fix the regular effective 
range via a truncation, we get the cut-off dependence 
found in Eq. ((95)) . Since L2{R) ^ R for regular solu- 
tions, we get an extra R^ suppression over the original 
one with one counterterm. Instead, if we use a Pade ap- 
proximant, Eq. (jl04p . we will get the convergence pattern 
of Eq. (|103[) . thus giving an extra suppression. 

The last case to consider is the one of an attractive 
singular potential. For it there is not a unique regular 
solution; the proper combination may be fixed by fixing 
the scattering length. Then the situation looks similar to 
the repulsive singular potentials. If we fix the effective 
range to the value obtained when fixing the scattering 
length, we will obtain an extra or R^ suppression 

depending on the usage or not of the Pade described in 



6. Overview of the Error Estimates and its Relation to 
Power Counting 

We can summarize the previous discussion about the 
truncation of the boundary condition in the following 
points 

• For an attractive singular potential it is necessary 
to fix at least one scattering observable in order to 
obtain a finite phase shift in the i? — )■ limit. 

• For a repulsive singular potential one cannot fix any 
scattering observable and at the same time obtain 
a finite phase shift in the i? — s- limit. This means 
that one is led to an alternative: either doing a 
finite cut-off computation or choosing the regular 
solution and remove the cut-off. 

• For a regular potential one can choose either to 
fix any scattering observable, and hence take the 
irregular solution, or not. 

• Fixing more scattering observables do not lead to 
a further (qualitative) short distance suppression, 
but a more singular potential (if attractive) im- 
proves the convergence. 

At this point we should stress the close connection be- 
tween the short distance suppression when one fixes a 
certain number of scattering observables and the issue of 
power counting in effective theories. For posing the dis- 
cussion in a more standard language, we will refer to the 
number of physical observables we are fixing as the num- 
ber of independent counterterms in the theory. It should 
be noted that fixing physical observables is not exactly 
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the same as setting counterterms. Counterterms are not 
observables, but they represent the unknown short dis- 
tance potential as a low energy expansion in terms of the 
delta function and its derivatives. On the other hand, 
setting physical observables is the visible effect for the 
need of counterterms. 

Although there are no clear rules in the literature, we 
will take the point of view that (i) the cut-off should 
be removed, (ii) it is the cut-off dependence what drives 
the construction of an acceptable power counting, and 
(iii) the long range potential is going to be treated non- 
perturbatively at any order (for a singular potential this 
is absolutely necessary). These assumptions have been 
adopted in our previous works 22, 25, 26] but it is fair 
to mention that they are not universally agreed upon, 
(i) is for example accepted in Ref. [1^, while rejected 
in Refs. [H, IH, [s^, (ii) is explicitly used in the RG 
analysis of Birse [H, [B, IH, and in Refs. [H [H for 
promoting certain counterterms, and (iii) is accepted 
in Refs. [481. J4p. [sot (the same which rejected (i)), but 
not in Ref. [23], in which it is advocated the perturba- 
tive treatment of the potential beyond LO, while LO re- 
mains non-perturbative (although no actual computation 
is done for NLO or NNLO). The amazing aspect of all 
these disagreements is that there is no operational def- 
inition of what would be a valid criterion and discrep- 
ancies are utterly based on the favourite prejudices of 
different authors (including ourselves). Here we are not 
going to discuss the validity or convenience of these as- 
sumptions (this has been partially done in our previous 
work [13, [2^ H^) but looking instead for their conse- 
quences. 

The first consequence of our previous analysis is that 
the long range potential and the counterterms (or short 
range potential) are not independent: the singularity 
structure of the potential determines whether countert- 
erms should be included or not. Only in the very special 
case of long range regul ar p otentials are the counterterms 
independent (see also [U[2^, [2^). This disagrees with 
power counting schemes based on naive dimensional anal- 
ysis, such as Weinberg's, in which all the counterterms 
with appropriate dimensions are included in the com- 
putation, regardless on the structure of the long range 
potential. But if the cut-off is to be removed, one must 
always include a counterterm if the long range poten- 
tial is singular attractive (in agreement with the previ- 
ous findings of Ref. [s^ and the results of Ref. [1^) and 
one cannot include any counterterm at all in case it is 
singular repulsive 



Actually, in momentum space one may find redundant countert- 
erms appearing in the short distance potential at fourth order 
in momentum Vsik',k) = Co + C2{k^ + fc'^) + C4,{k"^ + k'^) + 
Dik^k'^ -|- . . . . Obviously, C4 and -D4 are redundant or else 
do not correspond to the same order. This actually shows that 
there are may appear more counterterms than renormalization 
conditions. 



From a naive viewpoint, it might seem counterintu- 
itive that for a repulsive singular potential no countert- 
erm can be included, but as we discuss now it is indeed 
quite natural. If we consider a singular repulsive poten- 
tial with a characteristic long distance scale a, then the 
physics associated with a short distance scale as, such 
that as <C a, would not affect at all the long distance 
physics, since the long distance potential itself would act 
as a potential barrier which destroys any effect coming 
from scales 05 smaller than a. Then, there must be a 
strong short distance insensitivity which manifests itself 
as the dominance of the regular solution and thus a lack 
of counterterms when the cut-off is removed. 

The issue of repulsive singular potentials raises a un- 
expected consequence for the power counting of the long 
distance potential: if the potential is singular attractive 
at a given order, it should remain singular attractive at 
higher orders. If the long distance potential between two 
particles is known to be attractive, the effective theory 
should reproduce this feature to all orders in the expan- 
sion of the potential (if this is going to be used in any 
non-perturbative computation). Of course the full in- 
teraction between two particles cannot truly be singular 
attractive at all distances; if this were to be the case the 
system would collapse. However, in the effective theory 
it is the unknown short distance physics which are repul- 
sive, and these interactions are not explicitly modelled in 
effective theories, but implicitly via counterterms. 

Another important issue which arises from the previous 
RG analysis regards the number of counterterms and/or 
renormalization conditions one should include when one 
has an attractive singular potential. As we have seen, the 
qualitatively power law cut-off scaling behaviour does not 
depend at all on the number of counterterms included in 
the computation, but only on the power law divergence 
of the potential near the origin. Although this seems to 
suggest that there is no reason for adding counterterms 
beyond the first one, when we look at the quantitative 
cut-off dependence we can see that the situation may 
change. One clear example is given by the TPE potential 
at NNLO in the singlet channel. If we look at how 
much the phase shifts change between the cut-offs R = 
0.15 fm and R — 1.5 fm at a center of mass momentum 
k — 1.0 fm^^, we see 

\AS{k)\^(. ~ 25.21° (107) 
\A6ik)\^c ^ 15.11° (108) 

where the subscripts IC and 2C refer to the number of 
counterterms used. As one can appreciate, there is a no- 
ticeable improvement of the convergence when a second 
counterterm is added, and which can justify its inclu- 
sion It should be noted that this improvement is not 



As a matter of fact, it is interesting to notice that a third coun- 
terterm (to fix V2 = — 0.48fm) does not improve at all the con- 
vergence, but worsens it {\AS{k)\^i^ ~ 56.50°). 
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solely related with the second counterterm, but also with 
the fact that the effective range prediction when one fixes 
the scattering length, which is ro(ao) = 2.86 fm, is very 
close to the experimental value, rg — 2.77 fm, to which 
we fix the second counterterm. Taking into account that 
the convergence when fixing two counterterms depends 
on the A rg needed to correctly fix the effective range (if 
it were zero we would have an convergence pattern 
instead i?^) it is not a surprise that there is dramatic 
improvement. 

On the other hand, there are also good reasons not 
to include these extra counterterms. As argued in 
Refs. m, [2^ [2I], the inclusion of more counterterms 
than the minimum required for finiteness break orthog- 
onality constraints between different energy solutions, 
and in some cases, as the singlet at NLO also the 
Wigner causality bound [5l|, although the unphysical 
consequences of breaking this bound have not been stud- 
ied. Another good reason is the lack of an obvious im- 
provement in the NNLO results in the singlet channel, 
as can be seen in Fig. ([7|). While NLO phases notice- 
ably improve when fixing the effective range, the im- 
provement for NNLO phases is very small. Thus it may 
be a better option to preserve orthogonality. Actually, 
Fig. ([7]) suggests that it would be more instructive and 
perhaps profitable to improve on the long distance poten- 
tial than adding more and more counterterms Work 
along these lines is on the way [13] ■ 

It is also important to notice that the renormalizability 
of singular potentials depends on a very specific represen- 
tation of the short range physics, a Fade approximant if 
we write them in terms of an energy dependent bound- 
ary condition. This casts doubts on the renormalizabil- 
ity of momentum space treatments in which the contact 
interactions are parameterized as an ad-hoc expansion 
of deltas and their derivatives, since only a very precise 
short distance interaction will be able to renormalize the 
Lippmann-Schwinger equation. 



IV. BOUNDARY CONDITION AND VARIABLE 
PHASE EQUATION 

A. Equivalence between BC and Variable Phase 

The boundary condition, Eq. ([9]), corresponds to the 
inner radius i? of a boundary value problem defined in the 
region R < r < 00. We can give a physically appealing 
and computationally convenient interpretation of this BC 



The NLO prediction for the effective range is 2.29 fm [25| . smaller 
than the experimental one. 

This remark applies equally when going from a pure effective 
range expansion, which describes the low energy data with arbi- 
trary accuracy but randomly agrees to the intermediate energy 
region sensitive to the explicit pion exchange potential. 



in terms of a complementary outer boundary value prob- 
lem in the region < r < R. If we consider the family 
of potentials U{r, R) — U{r)6{R — r), which corresponds 
to a set of truncated potentials, ?7(r), at distances below 
a certain radius, r < R, acting only from the origin to 
the boundary radius R we would have at the boundary, 
r = R, the asymptotic wave function 

u{r)^sin{kr + 5{k,R)) r>R, (109) 

where we keep explicitly the dependence on R of the 
phase shift. The logarithmic derivative at the boundary 
r = R from the left is therefore 

LkiR) =kcot{kR + 6{k,R)). (110) 

If we identify this expression with that of Eq. ^ we get 
the following equation for S{k,R) 

Sli^ = _l{7(i?)sin2(fci? + 5(fc,i?)). (Ill) 
ait k 

This is a variable phase equation of the type ana- 
lyzed in Ref. [13], describing the evolution of the phase- 
shift, S{k,R), corresponding to the truncated potential, 
U{r,R) — U{r)9{R — r). The standard derivation for 
a regular wave function at the origin, u(0) = 0, is well 
known. For the case of general solutions including also 
energy dependence see e.g. [l^, [T5j. 

Thus, the renormalization group equation for the 
boundary condition of an inner truncated potential 
U{r, R) — U{r)6{r — R) is solved by the logarith- 
mic derivative of an outer cut-off potential U (r, R) = 
U{r) 9{R — r) through the variable phase equation 
Eq. (jllOp . This result builds a one to one relation be- 
tween the evolution of effective boundary conditions of 
the outer problem and that of a variable phases of the 
inner problem, which is illustrated in fig. ([8]). Obviously, 
the physical phase shift can be obtained from the variable 
phase as an asymptotic limit 

d{k)^6{k, 00). (112) 

On the other hand, the boundary condition extrapolated 
to the origin is given by 

Lfe(0+) = fc cot5(fc,0+). (113) 

In the standard variable phase approach [l3[ one assumes 
6{k,0) = corresponding to the absence of zero range 
interaction, and hence Lk(0) — 00, i.e. a regularity con- 
dition at the origin, u{0) = 0. In the present context 
it makes sense to define the short range phase-shifts as 
the variable phase extrapolated to the origin, once the 
potential has been completely switched off 

6s(k) ^ 6(k,0+) = lim S(k,R). (114) 

In fact, it turns out that the way the former limit must be 
taken is a bit subtle, particularly in the case of singular 
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FIG. 7: Phase shifts for the ^Sq single channel for the NLO TPE potential (left panel) and the NNLO TPE potential (right 
panel). We show the renormalized phase shifts computed when fixing one (IC curve) or two (2C curve) parameters, i.e. the 
scattering length and the effective range (ao = — 23.74 fm and ro — 2.77 fm). 
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FIG. 8: The relation between the boundary condition prob- 
lem defining the outer problem and the variable phase equa- 
tion defining the inner problem. The logarithmic deriva- 
tive at the boundary provides the identification Lk{R) ~ 
k cot{kR + 5{k,R)). 5{k,R) is the phase-shift produced by 
the potential(7(r) truncated at r > R. The physical phase- 
shift is given by 5{k) = 5{k, oo). The boundary condition at 
the origin is Lk{0'^) ~ kcotSs{k) where 5s{k) corresponds to 
the short distance phase shift, i.e., corresponding to U{r) — 
everywhere for r > (see main text). 



a long distance distortion of the short range phase shift 
Ssik) 

Ss{k) — > (Long range distortion) — > S{k) (115) 

In this interpretation the boundary condition regulariza- 
tion can be used to disentangle the short and long range 
physics from S{k) to obtain Ss{k) 

S{k) — > (Remove the distortion) — > Ss{k) (116) 

Another aspect of the variable-phase approach is the 
fact that we always deal with a given on-shell problem, 
that corresponding to the truncated potential. So, in 
the whole process there is no need to invoke any smooth 
off-shell behaviour, although we are changing the Hilbert 
space when the boundary radius is moved. 

For later purposes it is convenient to introduce the 
variable effective range M-matrix and its inverse the re- 
action V^-matrix, 



M(fc, R) 
V{k,R) 



k cot S{k, R) , 
tan 5{k, R) 



where direct insertion in Eq. Ijllip yields 



dM{k, R) 
dR 



UiR) 



and 

dV{k,R) 
dR 

respectively. 



= -U{R) 



M{k,R) 



sin kR 



sin kR 



cos kR 



V{k,R) coskR 



(117) 
(118) 

(119) 

(120) 



potentials (see also below). In the absence of a potential 
one gets a constant variable phase, and hence we would 
simply get S{k) — Ss{k). 

On these grounds, in the presence of a long range po- 
tential, the total phase shift 6{k) can be understood as 



B. Renormalization of Low energy Threshold 
Parameters 

Using the equivalence discussed above between the 
variable phase and boundary condition problems, it is 
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instructive to do a low energy expansion. At zero energy 
we have, 



ao{R) 



,. m 

am — ; — 

k^o k 



so one can obtain 



(121) 



(122) 



For increasing R this equation describes how the scat- 
tering length evolves as the potential U{r) is switched 
on. The physical scattering length is given by the value 
of ao{R) at R — oo, that is, ag = ao(cxD), provided we 
specify an initial condition at any point, say the origin. 
In the standard approach, one takes a(0) = 0. However, 
this is not the only possibility. We may use i? = cxd as 
the initial condition, and compute a(0) from there by 
integrating Eq. (|122p for decreasing R. This can be in- 
terpreted as the evolution of the scattering length as the 
potential is switched off for r > R. If we use an arbitrary 
value for a at i? = oo we obtain in general a{0) = as 0. 
Thus, we may interpret the value of as as the scatter- 
ing length corresponding to switching off the potential 
entirely. 

Going beyond the zero energy limit k — is not 
uniquely defined, because there are many equivalent ways 
of parameterizing the phase shifts by a low energy expan- 
sion. The coefficients of the expansion are, however, well 
defined. If for definiteness we use the effective range ex- 
pansion for the running Af -matrix 

M{k,R) = k cot S{k,R) 



ao[R) 2 



(123) 



It should be noted here that as, ro,S E^nd V2,s can only be 
defined for regular potentials The set of Eqs. (|124p . 
(|125p and (|126p for the running low energy threshold pa- 
rameters solve the set of equations for the parameters 
ioiR), and ^^{R), Eqs. ([71]), ([72]) and JISl) if in the 

relations ([501) and (|?T|) one substitutes the asymp- 

totic low energy threshold parameters for the running 
ones. 

The previous reasoning may also be applied to a low 
energy expansion of the variable reaction V^-matrix 



V{k,R) = 



ta.nd{k, R) 



^ao{R) - ^l3oiR)k^ + ■ ■ ■ (131) 



which generates the following equations for ao and (3q 

(132) 



^ = UiR){a,^R)\ 



dR 



-U{R) (R - ao) {R^ - 3aoR^ + 3/3o) (133) 



Obviously, we have the relation 

f3o{R) ^ roiR)aoiRf . 



(134) 



one has the set of equations [ij, [if 



It is straightforward to check that the set of Eqs. (|124p . 
and p25p and Eqs. (|132p and (|133p are mutually com- 
patible, as it should be. 

All low energy expansions share a common hierarchy 
for the low energy parameters; the evolution of a given 
low energy parameter contributing to a given order de- 
pends only on the evolution of lower order low energy 
parameters. The set of equations express the evolution 
of the low energy parameters at zero range when the long 
distance contribution is switched on. Conversely, they of- 
fer a possibility to determine the zero range low energy 
parameters from the total ones by downwards evolution 
in the cut-off variable R. 



2 1 



R_ 

ao 



(124) 



(126) 
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These equations have to be supplemented with the initial 
conditions [l^. [Tsj 

a(0+) = as a(oo) = a, (127) 

ro(0+) = ro,5 ro(oo) -ro, (128) 

1-2(0+) = Vo,S V2{(Xi)^V2. (129) 

(130) 



1. Application to NN Scattering 



The previous equations, Eqs. (fT24| . (fT25| and (fT26| . 
can be used to study the renormalization behaviour of 
the low energy parameters for neutron-proton scattering 
in the ^5o channel with the OPE potential. Although the 
OPE potential is regular for the ^5*0 channel, it behaves 
as 1 /r for distances below the pion Compton wave length. 
The consequence of this mild singularity is that the short 
distance scattering length, as = ao(0+) = 0, is zero, 
regardless the fact that we are not taking the regular 
solution of the OPE potential. This can be understood 



For attractive singular potentials ao{R), ro(-R) and V2{R) behave 
as highly oscillating functions for i? — > 0, so no definite limit can 
be obtained for R = 0. 
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by studying the behaviour of ao(-R) for the ultraviolet 
limit, R ^ (see Ref. [11]). We can solve Eq. pM]) 
for short distances and large scattering lengths with the 
Yukawa potential, given by Eq. (j57p . yielding 



^-c.ilogi? + C, 



(135) 



where C is an integration constant, and which is valid for 
small distances (mi? ^ 1) and large scattering lengths 
{R < ao{Rj). This solution of Eq. (|124p is equivalent to 
take the irregular solution of the Schrodinger equation or 
the stable ultraviolet fixed point studied in Sect. IIII D 21 
On the other hand, if we assume the scattering length to 
be small compared to the cut-off scale, aa{R) ^ R, we 
obtain 



ao{R) 



'2a 



(136) 



which turns out to be equivalent to take the regular so- 
lution of the Schrodinger equation or the unstable ultra- 
violet fixed point for the Yukawa potential. So the differ- 
ence between the regular and irregular solution lies in the 
ultraviolet behaviour of ao{R): although they both ap- 
proach the same ultraviolet limit, ao(O^) ~ 0, the trend 
is dissimilar. 

We can integrate Eqs. (HH]), (fHS)) and (fHS)) for the 
other low energy parameters of NN scattering in the ^ So 
channel with OPE, by taking as initial conditions for 
R = oo the experimental scattering length and effective 
range, while for the shape parameter V2 we use the value 
obtained from the Nijmcgen II potential, yielding 



a(0+) = 
ro(0+) = 4.46 fm 
i;2(0+) = 1.24 fm^ 



a(oo) = -23.74 fm, (137) 
ro(oo) = 2.77fm, (138) 
U2(oo) = -0.48 fm^ (139) 



In Ref. [T4!| we can find a more detailed discussion al- 
though different results are obtained due to the low in- 
frared cut used {Roo — 10 fm) In Ref. [53j the re- 
sult ro,s = 4.0 fm is obtained by a renormalization anal- 
ysis in momentum space, and in Ref. [Toj they obtain 
fQ.s — 3.10 fm although an ultraviolet cut-off of half the 
rho mass A = mp/2 is employed to perform the calcula- 
tion; if we use R — n/2A — n/mp ^ 0.8 fm, equivalent to 
the previous cut-off, we obtain tq.s — 2.86 fm. 



C. Low energy expansion of the short distance 
interaction 

So far, all we have done has to do with relating long 
and short range physics along the trajectory defined by 
the long distance potential. If it was for that nothing 
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would be achieved. We propose to make a low energy 
expansion of the short range physics. According to our 
previous discussion for a truncated potential it makes 
sense to expand either the M-matrix or the V^-matrix in 
powers of fc. 



Ms = k cot Ss = 



UkiO+ 

1 1,2 

as 2 
tanSs Ufc(0+) 



vsk 



k 

-as 



<(0+) 



1 



-Jsk^ - isk" + 



(140) 



(141) 



respectively. To achieve consistency with the low en- 
ergy expansions, Eq. (|45p and up to some order 
of the full M-matrix and ^-matrix we have to compute 
the short distance low energy parameters by integrat- 
ing downwards the set of Eqs. (|124p . (|125p and (|126p 
or Eqs. (|124p and (|133p . This way we exactly repro- 
duce the low energy expansion up to a desired order, 
and generate all higher orders in energy according to the 
long distance part of the potential. In particular, we also 
generate at any level of truncation the remaining higher 
order parameters. Thus, in the LO approximation we 
consider a and U{r) as independent parameters, corre- 
sponding to keep one term in Eq. p40p , and hence setting 
ro(0"*") — V2{Q'^) = • • • = 0. One then obtains ro, f2, etc., 
from a and U{r), and as a consequence the phase shift 
6(k). In the NLO approximation one considers a and 
To as independent variables and predicts V2,V3 etc., and 
hence the phase shift 5{k) from the knowledge of a, vq 
and C/(r). 

Note that if the physical low energy parameters, say 
a, ro and so on, and the potential U are known, the 
result is unique Another important point is that if we 
have a singular potential we cannot start directly at the 
origin, but at a given small radius R. An advantage of our 
method that will become clear below is that we can make 
R much smaller than any other scale in the problem. 
Even if we do this numerically, this is an effective way 
of eliminating the regularization. Moreover, if we take 
the limit of large R, much larger than the range of the 
potential a , we are effectively having a constant variable 
phase given by the short distance theory. It thus makes 
sense to compare the full result including the potential 
with that of the low energy expansion. 



: 20 fm is used instead. 



In practice this situation may be too optimistic, since low energy 
parameters may be directly deduced from the data, which are 
analyzed with a given model. In Ref. [sil we provide a thorough 
determination of the low energy parameters for NN interaction 
in all partial waves for the high quality potentials of Ref. [ssll . 
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V. HIGHER PARTIAL WAVES AND COUPLED 
CHANNELS 

The generalization of the present ideas to coupled 
channels is in principle straightforward and runs paral- 
lel to what was done in Sect. IIIll with some modification 
which we outline in the following. The coupled channel 
Schrodinger equation for the relative motion reads 



u"(r) + 



12 

u(-) + 



u(r) = k^u(r) , (142) 



where U(r) is the coupled channel matrix potential, u(r) 
is the reduced matrix wave function in the initial and final 
state and 



12 = diag(;i(/i + l),...,;Ar(/Ar + l))^ 



(143) 



= Aia.g{2^ll{E^El),...,2^XN{E-EN)), 

(144) 

are the angular momentum and the CM momentum in 
the coupled channel space respectively. Ei is the thresh- 
old energy and the reduced mass in the i-th channel. 
We assume for u(r) the boundary condition, 



W{R)-U{R)n{R)^(), 



(145) 



where (i?) is a real hermitean matrix in coupled chan- 
nel space, which in our framework encodes the unknown 
physics at distances r below the boundary radius R. In 
addition, we assume the asymptotic normalization con- 
dition for scattering states 



u(r) -> Uin(r) - Uout(»")S , 



(146) 



with S the standard coupled channel unitary S-matrix, 
S'l'S = 1, from which the scattering amplitude f = 1 + 
2ik^/^Sk^/^ can be obtained. The out-going and in-going 
wave functions can be defined as 

Uin(r) = h.^-\r), (147) 

Uout(r) = hW(r), (148) 
where h'^^-' and h'"' are given by 

h(+)(r) = diag(/i+(fcir),...,/i+ (fc^r)), (149) 

h(-)(r) = diag(/i;-;(fcir),...,V„(fcjvr)), (150) 

with hf{x) the reduced Haenkel functions of order /, 
hf(x) = xH^^i^{x) {fv^ = e^"), which satisfy the 
Schrodinger's equation for a free particle. 

To determine how should the matrix boundary condi- 
tion depend on the boundary radius in order to achieve 
the same S-matrix, we proceed similarly as in Sect. IIIII 
Making an infinitesimal displacement of the radius, R — > 
R + AR, and taking into account the total derivative of 
the wave function with respect to the boundary radius 



5u(r, R) 
OR 



URir, R) , 



(151) 



then, the derivative of the boundary condition is given 

by 

u"{R,R) + u'j^{R,R) ~ L',,iR)u{R, R)~ 

Lk{R)u'{R,R)-U{R)uR{R,R) = 0.(152) 

Deriving also Schrodinger's equation with respect to the 
inner radius R 



- u'^(r, R) + U(r)ui^(r, R) - k^UR{r, R) 
and the asymptotic wave function, Eq. (fT46l) 

dS 



UR{r, R) 
u'(r, R) 
u'«(r,i?) 



-uout(r)-^^, 

UinW - "out (OS, 

dS 



dR ' 



(153) 

(154) 
(155) 
(156) 



Thus, using Lagrange's identity we get 

= u(r,i?)^Ufl(r,i?)"-u"(r,i?)tufl,(r,i?), 
= {u{r,R)^URir,Ry -u'{r,RyuRir,R))' . 

(157) 

Integrating between R and oo and using the bound- 
ary condition, Eq. (|145p and Eq. (|157p and the identity 
u|^,jUout as well as the Wronskian u|JJUQ^^ — Uj^^Uout = 2ik 
we finally get 



2iStk^ =u(i?,i?)t 
dR 



k2 



U(i?)- 
L'{R) + L{Rf u{R,R). 



(158) 

This equation tell us how the S-matrix changes as the 
inner boundary radius is changed. If we require the S- 
matrix not to be dependent on the particular choice of R 
we get 



L(i?)2 + L'(i?) =U(i?) 



(159) 



which is the coupled channel generalization of Eq. 
which likewise accounts for the coupled channel momen- 
tum k dependent evolution of the boundary condition. 
The evolution of the low energy parameters can be trans- 
lated into the corresponding evolution of the short dis- 
tance boundary condition as a function of the boundary 
radius. Defining the dimensionless quantity 

Hk(i?) = RLkiR) = RKiR)ukiR)-^ , (160) 

and using Eq. (|159p we get 

i? ^(i?) = Hk(l - Hk) + U(i?)i?2 -f l2 - k^i?^ , 

(161) 
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a result already derived in our previous work [15| by dif- 
ferent means. Assuming for simplicity the degenerate 
case, ki = k and expanding into powers of the common 
momentum k one gets, 

Hfc(i?) = Ho(i?) + (fci?)2H2(i?) + ... (162) 

so the RG flow at zero energy becomes 

i?^ = Ho(l-Ho)+U(i?)i?2 + i2. (163) 

A thorough study for the general multichannel case is 
beyond the scope of the present work. Nevertheless, 
for some special cases some general conclusions may 
be drawn. The most appealing regime has to do with 
the possible appearance of chaos. If that would be so, 
the whole renormalization program would not be im- 
plementable in practice, since an absolute knowledge 
of the short distance conditions would be needed. It 
is well known, that the minimal order for a dynami- 
cal autonomous system to develop chaotic solutions is 
three since for two independent variables the Poincare- 
Bendixon theorem [s^ guarantees integrability. Actu- 
ally, if we introduce the variable t = — log(i?/a), as a 
new variable with i? as a dependent variable, the sim- 
plest case for potential scatte ring would correspond to 
two-coupled channels. In Ref. [15] we have discussed the 
case corresponding to OPE NN potential in the ^Si —^Di 
coupled channel, and infrared fixed points have been de- 
termined. 

To analyze the short distance behaviour, let us assume 
the potential to behave as an inverse singular power 



straightforward to obtain the equation for the variable 
S-matrix, 



U(i?)i?2 



a" 



(164) 



where a is a matrix with length dimension. We can diag- 
onalize the potential by a global transformation, say G, 
so that if we have a set of diagonal GCqG"^ at some 
initial value the solution will always be diagonal. This 
means that the system is integrable and hence chaos is 
precluded. 

Similarly to the findings of Sect. II VI fsee also Ref. fis'l) 
the coupled channel boundary condition, Eq. (|152p . 
for the outer boundary values problem, Eq. (jl42p and 
Eq. (|146p . can be interpreted in simple physical terms 
of a complementary inner problem where the potential 
U(r) acts in the interval R < r < oo. If we switch off 
the potential above a given boundary radius R we have, 
at the boundary 



Lfc(i?) - u'{R)u-\R) 

- K(i?) - <,,(i?)S(i?)] 

X [Ui„(i?)-Uo,t(i?)S(i?)]" 



(165) 



where S{R) is the S-matrix associated to the potential 
U(r) acting in the region < r < R, which inherits 
the dependence on the chosen boundary radius R. It is 



2ifc 



dS{R) 
dR 



S{R)h<-+HR)-h'-~\R) U(i?) 
h'^-\R)~h'^+\R)S{R)\ . (166) 



Further consequences of this equation, in particular its 
low energy limit can be looked up in Ref. [5J|. 



VI. SUMMARY AND CONCLUSIONS 

In the present work we have analyzed the role of bound- 
ary conditions in potential two-body scattering from a 
renormalization point of view. We remind that in quan- 
tum mechanical problems the most general form of in- 
teractions can be accommodated not only in terms of a 
given potential, but also as a nontrivial mixed bound- 
ary condition. The discussion has been carried out in 
configuration space, because there the disentanglement 
of the short and long distance physics is naturally for- 
mulated using both a suitable boundary condition at the 
origin and a local potential respectively. The suitability 
of the boundary condition is governed by a simple renor- 
malization group equation, and depends explicitly on the 
choice of the potential, the scattering energy as well as 
some renormalization conditions. In potential scattering 
the most appropriate conditions at low energies are the 
threshold parameters, such as the scattering length, the 
effective range and so on. The generalization of such an 
equation for the case of coupled channels is straightfor- 
ward and has also been analyzed in some detail. 

The resulting equations can be applied to many cases 
of interest, like the study of both infrared and ultravi- 
olet fixed points. We find stable and unstable infrared 
fixed points corresponding to the limit of small and large 
scattering lengths respectively, in agreement with previ- 
ous authors [3, • The same kind of fixed points is also 
found for regular potentials and repulsive singular poten- 
tials in the ultraviolet regime. In contrast, for the case of 
attractive singular power law potentials, we also describe 
ultraviolet limit cycles and attractors with a computable 
fractal dimension. Actually, these exotic renormalization 
group solutions are genuinely non-perturbative short dis- 
tance phenomena, and have a close relationship with the 
cut-off dependence of observables and hence with power 
counting. Once the RG trajectory of the boundary condi- 
tion enters the attractor, renormalized observables show 
a very smooth dependence with resp ect to the cut-off, 
as can be seen in Refs. [l^, [2^, l26j |. Moreover, these 
results imply fine tuning short-distance conditions as al- 
ready found in our previous works [3, fisj . Obviously, 
the appearance of an RG trajectory attractor and frac- 
tality immediately suggests the search for possible short 
distance chaotic behaviour of the renormalization group 
equations. Let us note that the presence of chaos at short 
distances would demand not only fine tuning but abso- 
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lute tuning of renormalization constants in the ultraviolet 
limit. Although our equations look very similar to those 
appearing in dynamical systems with chaotic behaviour, 
for the obvious candidate of coupled channel zero energy 
s-wave scattering with a singular power law potential we 
have not seen any traces of a chaotic pattern. 

Furthermore, the ultraviolet renormalization group be- 
haviour of singular potentials provides interesting in- 
sights into the power counting rules for non-perturbative 
effective field theories. Boundary conditions allow a clear 
estimate of the cut-off behaviour of scattering observables 
when certain renormalization conditions are imposed. 
Since one expects that an effective theory leads to a low 
energy description which becomes increasingly insensi- 
tive to detailed short distance dynamics, then the cut-off 
dependence of the observables can be used to construct a 
power counting for the unknown short range physics, i.e. 
the counterterms. One consequence is that counterterms 
are no longer completely independent of the long range 
potential, as happened in the original Weinberg's count- 
ing. Thus, a counterterm must be included in channels 
with an attractive sing ular potential, in agreement with 
previous works [H, [25l [39| . The inclusion of extra coun- 
terterms in these channels do not change the power law 
suppression in the cut-off of short range physics. How- 
ever, the proportionality coefficient does depend on the 
renormalization conditions. 

On the contrary in channels with a repulsive singular 
potential no counterterm can be added if the cut-off is 
to be removed [25j . This result, which seems very coun- 
terintuitive from a naive point of view, is indeed very 
natural when one takes into account that a long range 
repulsive potential acts as a potential barrier for the un- 
known short distance physics, thus completely screening 
their contribution to low energy physics. 

A very appealing perspective of the renormalization 
group method is provided by the definition of a comple- 
mentary scattering problem, which can be mapped into a 
variable phase equation [13] with nontrivial initial (short- 
distance) conditions and used extensively in our previous 
works [3, [H, [13] ■ Actually, using the long distance fixed 
points we find a kind of short-distance modified effec- 
tive range expansion both for large and small scattering 
lengths, which in some cases is amenable to a perturba- 
tive discussion. Moreover, our method allows to handle 
the case of singular potentials at the origin. This comple- 
ments the long distance modified effective range expan- 
sion proposed long-ago [31 and also deduced more re- 
cently by renormalization group arguments in Ref. [l^ . 
where by construction genuinely singular potentials at 
the origin where excluded. 

One of the most rewarding aspects of the previous 
and present investigations has to do with the practical 
elimination of short distance cut-offs in the NN scat- 
tering problem which long range pion exchange contri- 
butions become singular at the origin. Since Effective 
Field Theory ideas were proposed to study the NN prob- 
lem, a lot a progress has been made, but the existence 



of finite cut-offs has clouded the key renormalization is- 
sues. The renormalization program becomes necessary to 
make truly model independent calculations based on both 
the long distance physics explicitly governed by pion ex- 
changes and hence sensitive to chiral symmetry, and the 
unknown short distance dynamics. The present approach 
suggests a way from a renormalization group viewpoint 
to effectively remove these cut-offs, and suggests that the 
non-perturbative counterterms to achieve scale indepen- 
dence can indeed have a quite unexpected behaviour, as 
we have discussed and illustrated in detail for the ^So 
channel with the singular TPE Potential. Actually, our 
results are more general and can be applied in other con- 
texts where the unknown short distance physics plays a 
significant role. 
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APPENDIX A: SOLUTIONS TO THE 
SCHRODINGER EQUATION FOR SINGULAR 
POTENTIALS 

Singular Potentials are commonplace on Effective The- 
ories. Since the effective interaction is developed as a 
long range (or low energy) expansion, each new order 
gives rise to potentials which are more singular at short 
range scales. This can be seen in the NN chiral poten- 
tial, which at Q"^ order displays a l/r^'^'^ singularity for 
distances below the pion Compton wave length Con- 
sequently the study of the solutions of the Schrodinger 
equation for singular potentials is of main interest for 
this work, in which we treat the renormalization of the 
effective interactions in coordinate space. 

We will consider the case of a potential which shows 
a power law divergence near the origin, i.e. U{r) — > 



The degree of singularity can be milder in some cases, as in the 
^So singlet channel at order, in which the potential behaves 
as 1/r. 
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zba" ^/r". For it, we can write the Schrodinger equa- 
tion as 



u"{r) 



U{r) 



u = k u . 



(Al) 



For short enough distances the singular potential domi- 
nates and the centrifugal barrier becomes irrelevant; the 
conditions of applicability of the WKB approximation 
hold for r <C a (7i/2)^^*-"^^\ so we can use it to estimate 
the behaviour of the wave functions near the origin. For 
an attractive singular potential we have 



u{r) C y—j sin 



1 



n/2 - 1 



2/2-1 



■ V 



(A2) 



where ip is the semiclassical phase, which is a free param- 
eter that must be fixed from some given long distance 
information. It should be noted that there is no unique 
regular solution at the origin, as happened with regular 
potentials; both linearly independent solutions are regu- 
lar at the origin. On the contrary, for a repulsive singular 
potential we have 



u{r) 



' j'\ n/A 

C+[ - ] exp 
.a/ 



1 



n/2 - 1 Vr 



(D 



i/4 



exp 



n/2 



1 /a 



2/2-1' 



2/2-1' 



,(A3) 



which consists in a linear combination of an irregular (the 
positive sign exponential) and regular solution (the nega- 
tive sign exp.). Physical solutions demand that C+ = 0, 
thus precluding the possibility of fixing any observable 
via counterterms (see end of Sect. IIIIFj) . 

It is also interesting to consider the momentum expan- 
sion of the wave function 



Uk{r) = Uo(r) -f U2{r) + Ui{r) + 



(A4) 



which is used all along this work, specially when consid- 
ering the momentum expansion of the boundary condi- 
tion. The terms in these expansions obey the following 
differential equations 



U{r) 
U{r) 
U{r) 



l{l + l) 



1(1 + 1) 



"0 = 0, 



1(1 + 1) 



U2 



U4 



uo , 



U2 



(A5) 
(A6) 
(A7) 



and so on. Once a certain solution for uq (r) is given, the 
higher order terms can be solved via Green functions 



U2(r) 



dr'Ga(r,r')uo(r'), (A8) 



Ui(r) = / dr'Go(r,r')u2(r'), (A9) 
Jo 



where the Green function can be defined in terms of two 
linear independent solutions Uq^q and uq i, of the zero en- 
ergy Schrodinger equation 

Go(r,r') = Uo,a(r)uoMr')0(r' 'r) + 

uo,a(r')uQ^b{r)9{r~r'), (AlO) 



where uo.a and uq^i, fulfill the relation 
Uo,a(r)uoAr) " Uo,a(r-)uo,b W 



(All) 



corresponding to unity Wronskian normalization. From 
this representation of U2, u^, and so on, in terms of Green 
function it is very easy to obtain their power law be- 
haviour near the origin. Since any solution of the zero 
energy Schrodinger equation for a singular potential be- 
haves as r"/'', times some trigonometric or exponential 
function depending whether the potential is attractive or 
repulsive, we see that U2 behaves as r"/^+^ 



2/4 



U4 



„n+2 



as r""^^ X r"/'*, and so on, each power of k^ adding an 
^n/2+1 ^j^p wave function. This means that there is a 
strong short distance suppression in the expansion of the 
wave function in powers of momenta. 



APPENDIX B: EVALUATION OF THE 
TRUNCATION OF THE BOUNDARY 
CONDITION FOR ATTRACTIVE SINGULAR 
POTENTIALS 

In this appendix we want to discuss in further detail 
the error estimates when truncations are made in the 
boundary condition. Sect. IIIIFl thus breaking exact RG 
invariance. The idea behind these truncations is to fix 
some physical observables of the system, and then make 
predictions for the other observables. This violates exact 
RG invariance, which would imply a complete knowledge 
of the phase shifts of the system at all energies, some- 
thing which is not possible in practical computations, and 
which would also lack predictive power. Nevertheless, as 
we saw m Sect. HiTfI although exact RG invariance is 
broken for an arbitrary cut-off R, under some circum- 
stances we recover it in the i? — > limit. 

In all these cases, higher order RG equations are not 
fulfilled. Thus, if we only fix ao for example, we can 
rewrite the expansion of the boundary condition as 



Lk{R) = Lo{R) + k^ A2(R) 



(Bl) 



meaning by this notation that Lq fulfills its RG equation, 
while A2 and higher order terms don't. If we fix both ao 
and tq, we could write accordingly 



Lk(R) = Lo(R) + k^ L2(R) + Ai(R) 



(B2) 



in which now both Lq and L2 fulfill their respective RG 
equations, while the higher orders don't. In many occa- 
sion, the higher order terms can be neglected, and then 
the boundary condition reduces to 



Lk(R) 
Lk(R) 



Lo(R) 
Lo(R) 



fc2 L2(R), 



(B3) 
(B4) 
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which in turn impHes that A2/L0 or A4/L2 — > 
when the cut-ofF R goes to zero, thus simplifying the 
computations. 

For computing Lq, L2, and higher order terms in the 
expansion of the boundary condition it is useful to define 
the next momentum expansion of the wave function 



= uoir,R) + k^U2{r,R) + 
Ui{r, R) + ... 



(B5) 



in which we have normalized the wave functions accord- 
ing to the asymptotic normalization 



Uk{r,R) sin {k r + 6{k,R)) 



(B6) 



Uo{r,R) - 


r-ao, 


(B7) 


U2{r,R) - 


6 


(B8) 


Ui{r,R) - 







10 uq r^r — 120 W2 



(B9) 



and so on. Note that here S{k,R) is the phase shift de- 
pendence on the short distance cut-off (not to be con- 
fused with the variable phase). This normalization for 
the wave functions is consistent with the one used in the 
derivation of Eq. The terms uq, U2, etc, obey the 

following equations 



-u'^ + U{r)uo = 0, 



U (r) U2 
U{r) Ui 



U2 , 



(BIO) 
(Bll) 
(B12) 



which can be trivially deduced by applying the reduced 
Schrodinger equation to the expansion in powers of mo- 
menta of the wave function, Eq. (jBSp . 

Now let's study the case in which we fix the scatter- 
ing length of the system, thus fulfilling Eq. ((7T|) . while 
Eq. (|7^ and higher order ones are not. In such a case 
we can write the expansion of the boundary condition as 



LkiR) = LoiR) + k^ A2(i?) 



(B13) 



in which the relation of the terms in this expansion and 
the ones of the expansion of the wave function is given 

by 



LoiR) 

A2(i?) 



Uq 

W2U0 



U2Uq 



(B14) 
(B15) 



From the behaviour of the zero energy wave function for 
power-law singular potential, it is trivial to check that 
for small cut-offs 



LoiR) 



1 



cot 



2 ^r' 



(B16) 



On the other side, the behaviour of A2 can be easily 
evaluated if we take into account the following Lagrange 
identity 



(M2U0 



U2W0)' 



(B17) 



where we have dropped the dependence of uq on the cut- 
off, since this wave function is cut-off independent. From 
this we find that 



A2(i?) 



1 



uoiR)- JO 

which for small cut-offs approximately behaves as 

2 



A2iR)^Rf 



2 ^r' 



fo 



(B18) 



(B19) 



with f{x) some unspecified trigonometric function. It 
can be easily shown that A2 obeys the differential equa- 
tion 

Ai(i?) + 2Lo(i?)A2(i?)-l = 0, (B20) 

in contrast with the rcnormalization group equation 

L'2{R) + 2 Lo{R) L2{R) + 1 ^ . (B21) 

As we can see, fc^ A2 <C Lo, which means that we can 
drop the A2 term, without jeopardizing the smooth cut- 
off dependence of the phase shifts for small cut-offs. In 
fact, it is curious to see that dropping A2 improves the 
convergence. When we take Lk{R) = Lo{R), we arrive 
at the following cut-off dependence for the phase shift 



^=-k\ 
dR V 



fuk{R, R) 



(B22) 



while when considering the contribution for the A2 term, 
i.e. we take Lk{R) = Lo{R) + k^ A2{R), we arrive at the 
following result 



dS 
dR 



-2t' 



UkiR,R) 



k^R"-^"^ , (B23) 



which is qualitatively the same dependence, although the 
cut-off dependence is doubled. In short, what we can 
see is that, when fixing the scattering length, A2(i?) is 
negligible in comparison with Lq{R), and can therefore 
be removed. 

Unfortunately this changes when considering the the- 
ory in which one fixes both the scattering length and 
the effective range. In it, the boundary condition can be 
expanded as 



Lk{R) = Lo{R) + fc2 L2{R) + k^ Ai{R) 



(B24) 



and, as we will show below, the A4 contribution cannot 
be ignored. The behaviour of Lq{R) is the one given 
by Eq. (|B16p . i.e. exactly the same as in the previous 
case. On the contrary, the behaviour of L2{R) differs 
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significantly from that of A2(i?), given by Eq. (|B19[) . 
The reason hes in the fact that we are fixing the effective 
range. This can be deduced from the behaviour of the U2 
wave function, which obeys the equation 



1*2 + U (r) U2 Uq . 



(B25) 



The solutions of the above equation can be separated into 
an homogeneous and inhomogeneous piece 



U2 = U2,H + U2,I , 



(B26) 



which in turn means that L2{R) can be written as the 
contribution of these two pieces 



L2 = L 



2.H 



where 



L2J 



(B27) 

(B28) 
(B29) 



The inhomogeneous piece of U2 comes from the nontrivial 
contribution stemming from the uq wave function and 
can be computed via Green's functions. Its short distance 
behaviour is much smoother than that of uq, U2j scales 
as r^"/''+^ times an oscillating function f{x) with x — 
("h)"^^"^' ^hil^ ""o just scales as r"/*. From these 
behaviours it is trivial to see that 



L2j{R)^A2{R)^R, 



(B30) 



in which we wanted to make clear the identification be- 
tween the inhomogeneous piece of L2 and what we pre- 
viously called A2. 

The problem with the inhomogeneous piece of U2 is 
that its contribution to the effective range is fixed, and 
is given by the following formula 



^0 "'0 



dr [{r - aof - ulir)] , (B31) 



as can be trivially deduced from Eq. (|B17p . and in which 
we have taken the cut-off to zero for simplicity. For mod- 
ifying the total effective range we need to add one contri- 
bution from the asymptotic behaviour at large distances 
of the homogeneous wave function in such a way that 



U2,H 



■ ao ro^H r . 



(B32) 



Then, the complete effective range of the wave function 
1*2 = 1*2, _ff + 1*2, / will be Tq = H + tq j, thus fixing 
the effective range to the desired value. But the homo- 
geneous contribution to U2 is much more singular than 
the inhomogeneous one. The short distance behaviour of 
U2,H is given by 



U2,HiR) R"^'^ sin 



(B33) 



with (f2 a semiclassical phase. From the previous be- 
haviour for 1*2, ij, Eq. (|B33[) . one can see that L2.H (and 
therefore L2) will behave in a very similar way to Lqj 
i.e. a singularity times some oscillating function. 

But one can get a much better evaluation of L2^h by 
considering the following Lagrange identity for uq and 

U2.H 



("2,_ff"0 - W2,Hlio)' = > 



(B34) 



in contrast with U2j which obeys the Lagrange identity 
Eq. (|B17p . This means that u'2 j^uq — i*2,_ffi*o, which ap- 
pears in the numerator of L2^h, is a constant value. Eval- 
uating at large distances, we obtain 



'U2,hUq - U2,HUo 



from which trivially follows 



fo,H 2 



L2,h{R) 



"0 \ ra^H 



uoiR) 



(B35) 



(B36) 



It is curious to see the equations which both contribu- 
tions to L2 follow 

4,^ + 2Lo(i?)i2,H(i?) + 2 - 0, (B37) 
L'2j + 2Lo{R)L2j{R)-l = 0, (B38) 

from which taking into account that L2 = L2.H + L2J 
one recovers the renormalization group equation. 

Finally, to evaluate A4 we need its expression in terms 
of wave functions 

A4(i?)^ "^"°"2"^"° +^L2(fl), (B39) 
1*0 Wo 

where the 1*4 contribution stems solely from the inhomo- 
geneous piece (since we are not fixing W2). By considering 
the appropriate Lagrange's identity 

(U4U0 - U4W0)' = "0 "2 , (B40) 

it is trivial to deduce the entire evaluation of A4(i?) 

1 

A4(i?) = ^/ druoU2 + —L2[R). (B41) 
"o Jo "0 

Due to the U2L2/UQ contribution, we see that A4 is as 
singular as L2, and therefore as Lq, so its presence cannot 
be neglected in the momentum expansion of the bound- 
ary condition. The same happens with Ag, Ag, and so 
on, but for the purposes of this appendix is enough to 
consider only A4. 

With all this we are prepared to evaluate the cut-off 
dependence of the phase shifts, which is given by 



dS 



— = -k' [Ll + A'^ + 2LoAi] ul{R) 



(B42) 



28 



After evaluating A4 inside the brackets, we are left with 
the much simpler expression 



dR 



-/c^ U2(i?) uo{R) + 0{k'^) 



(B43) 



which is just the same cut-off dependence as in the pre- 
vious case, but somewhat diminished by the appearance 
of an extra fc^ factor. 

Of course this is only the first step in showing that the 
phase shift is well behaved. The 0{k'^) term needs the 
explicit consideration of Ag, the 0{k^) term of Ag, and so 
on. Although it is straightforward to show that including 
the appropriate expansion of the boundary condition, the 
phase shift has a convergent behaviour for i? — > 0, it 
is very inconvenient for doing practical computations in 
which we want an specific expression for the boundary 
condition at the boundary radius R. 

The solution to these problems is relatively straightfor- 
ward if we realize that the contribution to A4, Ag, etc, 
which cannot be ignored comes from the homogeneous 
contributions to the wave functions, i.e. from uq and U2 
(when we fix ao and tq). For the general case in which 
we fix ao, To, 1)2, ■•• , I'n, the correct way to truncate the 
boundary condition is 



Lk{R) 



„2 



U2n , 



Uq + fc^ W2 + ■ ■ • + fc^" U2r, 



A , (B44) 



where A contains the inhomogeneous higher order pieces. 
It is easy to see that A(i?) ~ R, so it is negligible with 
respect to the main piece, which behaves as ~ 
The main piece generates all the nontrivial terms needed 
for having a well-behaved R ^ limit on the phase shift. 
In fact, this probes the complete uniqueness of the result 
when fixing, for example, either ao and tq in the k cot 6 
expansion, or ao and (3o in the tanS/k expansion. 



APPENDIX C: SOFT AND HARD BOUNDARY 
CONDITIONS 

In this appendix we analyze how can we define contact 
interactions through boundary conditions and the rela- 
tion between what we call soft and hard boundary con- 
ditions. Naively, the scattering amplitude corresponding 
to a contact interaction is given by 



fik) 



1 



— 1/a — ik 



(CI) 



This amplitude always has a pole at k — i/a which cor- 
responds to a virtual state for a < and to a bound 
state for a > The limit of trivial scattering a ^ 



There are situations where this pole is absent. For instance a 
square well repulsive potential 



corresponds to send the pole, either positive or negative, 
to infinity. Or rather, to take S = mr. By continuity of 
5{k) in a we should have a difference of vr if a — s- 0+ or 
if a — > 0~, thus the limit of trivial scattering a is 
not uniquely defined. We can obtain the amplitude (|C1[) 
by solving the free Schrodinger equation using the "soft" 
boundary condition at the origin 



au'(O) + m(0) = 0. 



(C2) 



If we take the limit a ^ we get the "hard" condition 
u{0) = 0. This gives a trivial phase-shift 6 = 0. Note 
that this limit is indeed singular since it does not distin- 
guish between a — > 0"*" and a — > 0". So, how can we 
have nontrivial scattering with no pole? As discussed in 
Ref. [11] there is no way of doing this with a finite range 
potential in the limit of vanishing range. 

The proper way is to use a "hard" condition at the 
point r = a for a > 0, 



yielding 



fik) 



u{a) = , 



1 



fccot(— afc) — ik 



< a < 00 . 



(C3) 



(C4) 



This amplitude does not have any pole in the complex 
plane. Thus, one starts from a = 00 to a = 0^ using 
Eq. (|C3[) which corresponds to a repulsive core and then 
goes on for a = 0~ with Eq. (jCSj) . 

For the purely short distance theory, the generaliza- 
tion to coupled channel scattering is almost trivial. If 
we assume s— wave scattering, and using the boundary 
condition at the origin 



au'(0) + u(0) = 0. 



(C5) 



with a the s-wave scattering length matrix, we get for 
the coupled channel amplitude 



ik. 



(C6) 



with k given by Eq. (|144|) . The simplicity of the 
derivation contrasts with the cumbersome treatment of 
the Lippmann-Schwinger equation with a sharp three- 
momentum cut-off presented in Ref. [5§|. 



APPENDIX D: WHY A BOUNDARY 
CONDITION INSTEAD OF A POTENTIAL AS A 
REGULATOR IN COORDINATE SPACE ? 

In this appendix we elaborate further on the suitabil- 
ity of the BC as compared to the use of square well or 
delta shell regulators. Our main concern has to do with 
choosing a regularization method where a truncated low 
energy expansion of the amplitude corresponds exactly 
with a truncated low energy expansion of the regulator. 
One of the advantages of a boundary condition as a short 
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distance regulator, as opposed to a short range potential, 
is related to the non appearance of spurious high order ef- 
fects in a low energy expansion of the amplitude . As we 
will see here this requires infinitely many counter-terms 
in an energy dependent potential. Of course, any of these 
methods when considered to all orders ought to yield the 
same result; the difference has to do with truncating the 
regularized potential to a given order. As we will show 
the BC is the only method of the three analyzed where 
truncation is consistent order by order. To understand 
the situation let us consider the three regularizations in 
coordinate space: Square well regularization |3^, delta 
shell re gula rization [l7| and boundary condition regular- 
ization [ij]. 



Square well regulator. In this case we get 

Uk{r) = {Uo + k^U2 + ...)eia-r). 
The phase-shift reads 



k cot S — 



VF + C/fc cot(fca) cot{y/W+Uka) + k 
kcot{ka) - Vfc^ + Uk cot(Vfc2 + Uku) 
Thus, the scattering length is given by 

tan(A/c7oa) ' 



/Uoa 



(Dl) 



(D2) 



(D3) 



and the effective range 



ro 



2a 1 



1 



alloa 



a 

3^ 



a-a / 



Uoa + aaUo-Uoa-"] (D4) 



and so on for V2, etc. Obviously, we can always fit 
any set of low energy parameters by adding suffi- 
cient energy dependent terms to the potential Uq, 
U2 etc. any value of the low energy constants we 
want. The mapping, however, is not one to one; we 
have infinitely many solutions for Uq Eq. (|D3p for 
a given value of a at a given scale a. This multi- 
valuation propagates to higher order low energy pa- 
rameters. As a consequence, if we only specify a set 
of low energy parameters, the remaining higher or- 
der ones are multivalued, and hence the phase-shift 
is not uniquely predicted. To make the ambiguity 
more explicit let us take the limit a ^ and con- 
sider the LO truncated potential, Uq. In this limit 
we obtain from inversion of Eq. (jD3p , 



(D5) 



with n an arbitrary integer number. Using this 
asymptotic solution we get for the effective range. 



ro = 2a 1 



1 



a(n-M/2)27r2 



(D6) 



which depends manifestly on the arbitrary value of 
n. The n dependence might be cancelled by in- 
cluding a NLO term, 1/2. More generally, for any 
finite a one should include infinitely many terms in 
the energy expansion of the square well potential, 
Eq. (|D1[) to get rid of these multi-valuation in the 
pure short range theory. 



• Delta Shell regulator. In this case we get 
U{r) = {Uo + k^U2 + ...)5{a-r). 
Then one gets for the scattering length 



(D7) 



(D8) 



(D9) 



1 + ?7oa2 ' 

whereas the effective range reads, 

4a / a \ I a ^ 

ro = -(l-— +2aC/2(l- 

3 V 2q/ V a/ 

In this case there are no multiple solutions, but 
similarly to the square well regularization we get 
spurious terms at higher orders 

Boundary condition regulator. In this case we have 



= fccot(fca -f (5) = $0 + ikaf^2 + 



<■(«) 
life (a) 

In this case the translation between the bound- 
ary condition and an effective range expansion is 
most straightforward; each energy contribution to 
the BC provides a new term in the effective range. 
If we take the limit a ^ we get a one-to-one 
mapping. Thus, the BC method provides a com- 
patible hierarchy of equations in a low momentum 
expansion. 



APPENDIX E: RENORMALIZATION GROUP 
EQUATION FOR POTENTIAL 
REGULARIZATION 

The standard way of visualizing renormalization is by 
means of potentials. In this section we use our renormal- 
ization group ideas similar to those presented in Sect. Hill 
to determine the RG evolution of these short distance 
regulators. 



1. Square well Potential Regularization 



Actually, one of the advantages of using Dimensional regulariza- 
tion (DR) in the MS scheme was that there was a one-to-one 
connection between the bare potential and the renormalized am- 
plitude. Unfortunately, DR has never been applied to the non- 
perturbative regularization of singular potentials at the origin. 



In Refs. [39j a short range energy dependent square 
well has been employed to regulate the short distance 
behaviour. That means taking the family of potentials 

Uk[r,R) - Uk{R)Q{R-r)^U(r)d(r~R). (El) 
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The regular solution for r < R is (assuming Uk{R) < where A has to be determined by matching to the long 
for definiteness) range piece of the wave function. Then 



u{r) = A sin (^^Jk^ - Uk{R)r^ r<R, (E2) 



ldS_ 

k dR 



sin2 ^k^-Uk{R)R) {Uk{R) - U{R)) + U'kiR) 



^ sin {2^k-^-Uk{R)R 
2" 4VA;2 - Uk{R) ^ 



(E3) 



If we demand the phase shift to be independent on 
the short-distance regulator we get the renormalization 
group equation, 



KU'^,{R) = {Uk{R)-U{R)) 



2sin2(,/)) 



1 — cos sin (j) 
where we have defined the dimensionless combination, 



(E4) 



^k^ -Uk{R)R. 



(E5) 



A similar equation has also been found in Ref. [4l| by 
different means. The low energy fixed points are given 
by 



sin(0) = 1 = 0sin(0) cos((/)) 



(E6) 



The first equation has the analytical solution 0„ = nir. 
The other equation has also infinitely many solutions 
All these fixed points are stable, so we have infinitely 
many branches, in agreement with the observation of 
Ref. fH. 



2. Delta shell Potential 

In Ref. [It] a delta shell potential regularization for 
the short range potential has been introduced, and a RG 
equation in momentum space has been obtained by cut- 
ting off the high energy components. According to our 



point of view, we should cut-off the long range potential 
at short distances. Then, we have 

U{r) = Uk{R)S{r R) + U{r)e{r - R) . (E7) 

Let us denote by and VLir) the wave functions regular 
and singular at the origin respectively associated to the 
long range potential alone, U{r), and fulfilling 



UL sin(fcr + Sl) , 
VL — > cos(fcr + Sl) ■ 



(E8) 
(E9) 



If we take the u(r) = AuL(r) + BvL{r) for r > R, then 
we have 



u'jR) 
u{R) 



kcot{kR) = UsiR) . 



(ElO) 



To obtain the RG equation let us compute the change of 
the phase shift. 



AS 



^-A [Us{R)u{R)^] + \uL{R)ARu{Rf .(Ell) 



k L -V ^ V ^ J ■ 

Using the condition lElOl we get 

U's{R) = Ul{R) - 2 [Us{R) + kcot{kR)] Us{R) .(E12) 
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